Scattering, determinants, hyperfunctions in relation to 



The method of realizing certain self-reciprocal transforms as (absolute) scattering, previ- 
ously presented in summarized form in the case of the Fourier cosine and sine transforms, is 
here applied to the self-reciprocal transform f(y) i— > Tt(f)(x) — J* °° Jo{2y/xy)f(y) dy, which 
is isometrically equivalent to the Hankel transform of order zero and is related to the func- 
tional equations of the Dedekind zeta functions of imaginary quadratic fields. This also allows 
to re-prove and to extend theorems of de Branges and V. Rovnyak regarding square inte- 
grable functions which are self-or-skew reciprocal under the Hankel transform of order zero. 
Related integral formulae involving various Bessel functions are all established internally to 
the method. Fredholm determinants of the kernel Jo(2^fxy) restricted to finite intervals (0, a) 
give the coefficients of first and second order differential equations whose associated scattering 
is (isometrically) the self-reciprocal transform Tt, closely related to the function -fej ■ Re- 
markable distributions involved in this analysis are seen to have most natural expressions as 
(difference of) boundary values (i.e. hyperfunctions.) The present work is completely indepen- 
dent from the previous study by the author on the same transform TC, which centered around 
the Klein-Gordon equation and relativistic causality. In an appendix, we make a simple-minded 
observation regarding the resolvent of the Dirichlet kernel as a Hilbert space reproducing kernel. 
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1 Introduction (the idea of co-Poisson) 

We explain the underlying framework and the general contours of this work. Throughout the paper, 
we have tried to formulate the theorems in such a form that one can, for most of them, read their 
statements without having studied the preceeding material in its entirety, so a sufficiently clear 
idea of the results and methods is easily accessible. Setting up here all notations and necessary 
preliminaries for stating the results would have taken up too much space. 

The Riemann zeta function £(s) = jj + ^ + gj + . . . is a meromorphic function in the entire 
complex plane with a simple pole at s = 1, residue 1. Its functional equation is usually written in 
one of the following two forms: 

^ f r(^)CW-^^ £ r(i T ^)C(i-.s) (ia) 



C(*)=Xo(*)C(l-«) *>(*)= Tr-*^^ (lb) 



r(ig) 
r(f) 



The former is related to the expression of tt i r(|)£(s) as a left Mellin transform 1 and to the Jacobi 
identity: 

r(-)C(s) = g J Q (0(t) l)**" 1 dt (»(*) > 1) (2a) 

1 f°° 1 

= -J o (6(t)-l--^)t*- 1 dt (0<B(*)<1) (2b) 

e(t) = i + 2j2q n2 i = e-" t e{t) = ±e{\) (2c) 

n>l 



1 in the left Mellin transform we use s — 1, in the right Mellin transform we use —s. 
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The latter form of the functional equation is related to the expression of £(s) as the right Mellin 
transform of a tempered distribution with support in [0, +00), which is self- reciprocal under the 
Fourier cosine transform: 2 



/>oo 

/ (£u 



COO = / (Y6 m (x)-l)x- s dx (3a) 



poo 

/ 2 cos(2nxy)C^2 S n(y) -l)dy= Y] (5 m (x) - 1 

JO 



(a > 0) (3b) 



n>\ m>l 

This last identity may be written in the more familiar form: 

/e 3 «^i„(#=^S m (i) (4) 

which expresses the invariance of the "Dirac comb" distribution ^2 m ^S m (x) under the Fourier 
transform. As a linear functional on Schwartz functions <fi , the invariance of ^2 m £%,8 m (x) under 
Fourier is expressed as the Poisson summation formula: 

£ 0(n) = *M = / e 2 ""^^) dx (5) 

The Jacobi identity is the special instance with <fi(x) = cxp(— irtx 2 ), and conversely the validity of 
(5) for Schwartz functions (and more) may be seen as a corollary to the Jacobi identity. 

The idea of co-Poisson [4] leads to another formulation of the functional equation as an identity 
involving functions. The co-Poisson identity ((10) below) appeared in the work of Duffin and 
Weinberger [13]. In one of the approaches to this identity, we start with a function g on the positive 
half-line such that both g(t)dt and L°° g(t)i _1 dt are absolutely convergent. Then we consider 
the averaged distribution g * D(x) = J °° g(t)D(j) * where D{x) = J2 n >i $n( x ) — l x> o(x). This 
gives (for x > 0): 

i.aw.fiWa.rsM* (6 ) 

n=l n J ° 

If g is smooth with support in [a, A], < a < A < 00, then the co-Poisson sum g * D has Schwartz 
decrease at +00 (easy from applying the Poisson formula to g ^Y ; cf. [8, 4.29] for a general 
statement). The right Mellin transform g * D(s) is related to the right Mellin transform g(s) of g 
via the identity: 

/•oo />oo 

<fTD(s) = (g*D) {x)x~ s dx = C(s) / (aOaT a = ((s)g( S ) (7) 
Jo Jo 

This is because the right Mellin transform of a multiplicative convolution is the product of the 

right Mellin transforms. The necessary calculus of tempered distributions needed for this and other 

statements in this paragraph is detailed in [8]. The functional equation in the form of (lb) gives: 3 

g^D( S ) = X o(s)C(l ~ s)g(s) = X o(s)l(gJ^D(l - s) I(g)(t) = 9 -^p- (8) 



2 of course, <5 m (x) = S(x — m). 

3 one observes that I(g)(s) = g(l — s). 
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One may reinterpret this in a manner involving the cosine transform C acting on L 2 (0, +00; dx). 
The Mellin transform of a function f(x) in L 2 (0,oo;dx) is a function f(s) on 0?(s) = i which is 
nothing else than the Plancherel Fourier transform of e^ u f(e u ): /(| + 17) = / °° f(x)x~^~ 11 dx = 
fZefWete-^du, J °°\f(x)\ 2 dx = ,Tj/(e")ei| 2 <fo = i / R(>)= x \f(s)\ 2 \ds\. The unitary 
operator CI is scale invariant hence it is diagonalized by the Mellin transform: CI(f)(s) = Xo( s )/( s ), 
C(/)( s ) = Xo(s)/(l — s )i where Xo(s) is obtained for example using f(x) = e~ 7TX and coincides with 
the chi-function defined in (lb). It has modulus 1 on the critical line as C is unitary. So (8) says 
that the co-Poisson intertwining identity holds: 

C(g* D) = 1(g)* D (9) 

The co-Poisson intertwining (9) or cxplicitely: 

r 2 cos(2 7 r^)ff;^M- r im dt \ dx= f^iMri_ r g{t)dt iy>0) (10) 

, ' \m=l m Jo t J n=1 U Jo 

is, when g is smooth with support m [a, A], < a < A < 00, an identity of (even) Schwartz 
functions. If g is only supposed to be such that J °° |g r (t)|(l + j) dt < 00 then the co-Poisson 
intertwining C(g * D) = 1(g) * D holds as an identity of distributions (either considered even or 
with support in [0,oo)). Sufficient conditions for pointwise validity have been established [8]. The 
general statement of the intertwining is C(g * E) = 1(g) * C(E) where E is an arbitrary tempered 
distribution with support on [0, 00) (see footnote 4 ) and it is proven directly. The co-Poisson identity 
(10) is another manner, not identical with the Poisson summation formula, to express the invariancc 
of D under the cosine transform, or the invariance of the Dirac comb under the Fourier transform. 

If the integrablc function g has its support in [a, A], < a < A < 00, then g * D is constant in 
(0, a) and its cosine transform (thanks to the co-Poisson intertwining) is constant in (O,^!" 1 ). Up 
to a rescaling we may take A = a" 1 , and then a < 1 (if a non zero example is wanted.) Appropriate 
modifications allow to construct non zero even Schwartz functions constant in (—a, a) and with 
Fourier transform again constant in (—a, a) where a > is arbitrary [8]. 

Schwartz functions are square-integrable so here we have made contact with the investigation of 
de Branges [1], V Rovnyak [28] and J. and V. Rovnyak [29, 30] of square integrablc functions on 
(0,oo) vanishing on (0,a) and with Hankel transform of order v vanishing on (0,a). For v = — ^ 
the Hankel transform of order v is f(y) 1— ► / °° cos(xy)f(y) dy and up to a scale change this 
is the cosine transform considered above. The co-Poisson idea allows to attach the zeta function 
to, among the spaces defined by de Branges [1], the spaces associated with the cosine transform: 
it has allowed the definition of some novel Hilbert spaces [3] of entire functions in relation with 
the Riemann zeta function and Dirichlct L-functions (the co-Poisson idea is in [4] on the adeles of 
an arbitrary algebraic number field K\ then, the study of the related Hilbert spaces was begun for 
K = Q. Further results were obtained in [7].) 

4 both sides in fact depend only on E(x) + E(—x) as a distribution on the line, which may be identically 0, and 
this happens exactly when E is a linear combination of odd derivatives of the delta function. 
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The study of the function xo{s) = n s ~2 r(i) ' °^ un ^ m °dulus on the critical line, is interesting. 
We proposed to realize the xo function as a "scattering matrix" . This is indeed possible and has been 
achieved in [6]. The distributions, functions, and differential equations involved are all related to, 
or expressed by, the Frcdholm determinants of the finite cosine transform, which in turn are related 
to the Fredholm determinants of the finite Dirichlct kernels sln ( t ( a: ~^)) on r_i \] The study of the 
Dirichlct kernels is a topic with a vast literature. A minor remark will be made in an appendix. 

i rf 1—5 ) 

We mentioned the Riemann zcta function and how it relates to xo{ s ) = 7T S ~2 r(l) an d to the 
cosine transform. Let us now briefly consider the Dedekind zeta function of the Gaussian number 
field Q(i) and how it relates to x( s ) = r[7) ano - to the H transform. The H transform is 

H(g)(y)= J (2^xy)g(x)dx j ( 2<v /5y) = (11) 

Up to the unitary transformation g(x) = (2x)~ ^ f (V2x) , H(g)(y) = (2y)~ik(^2y), it becomes 
the Hankel transform of order zero k{y) = ^/xyJo(xy)f(x) dx. It is a self-reciprocal, unitary, 
scale reversing operator (Tt(g(\x))(y) = jT~L{g)(j))- We shall also extend its action to tempered 
distributions on R with support in [0, +oo). At the level of right Mcllin transforms of elements of 
L 2 (0,oo;dx) it acts as: 

H(g)(s)=x(s)g(l-s) x ( fl ) = l£L=- f) R( a ) = I (12) 

T(s) 2 

It has e~ x l x >o{x) as one among its self-reciprocal functions, as is verified directly by series expansion 
f °° J (2^xy)e~y dy = £^ =0 { -^fx n y n e^ dy = e~*. The identity 

J (2Vt)t- s dt = X (s) = ^ (13) 

r(s) 

is equivalent to a special case of well-known formulas of Weber, Soninc and Schafhcitlin [33, 
13.24.(1)]. Here we have an absolutely convergent integral for | < 5R(s) < 1 and in that range 
the identity may be proven as in: e~ x = Jo(2y/xy)e~ y dy = J °° Jo(2y/y)-^e~i dy, T(l — s) = 
Jo°° "^o(2y / y)(/ ° x~ s ~ 1 e~* dx) dy = T(s) Jo(2y/y)y~ s dy. The integral is semi-convergent for 
5R(s) > 4, and of course (13) still holds. In particular on the critical line and writing t = e u , 
s = | + ij, we obtain the identities of tempered distributions / R e^ u Jo(2ei u )e~ nu du = x{\ + *7)i 
ei"J (2ei«) = £/ H x(| +il)e +i ^du. 

We have Co(i)(«) = \ E(n,m)#(o,o) ( n *+ m 'y = + £ + h + h + h + ' ' ' = E n >i ^ and !t 
is a meromorphic function in the entire complex plane with a simple pole at s = 1, residue ^. Its 
functional equation assumes at least two convenient well-known forms: 

(V4r(2^rr( S )c QW ( S ) = (v^-^rf^ra - S )c m (i *) (u a ) 

(^) S C Q(4 )( S ) = XW(^) 1 - S C QW (1 - s) X (s) = (14b) 
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The former is related to the expression of ir s T(s)(q^(s) as a left Mellin transform: 

1 1"°° 

^~ s r(s)CQW (s) = i J (0(t) 2 - l)* 8 " 1 dt (&(*) > 1) (15a) 

= - / (6(t) 2 - 1 - (0 < K(s) < 1) (15b) 

9(tf = U{\f (i5c) 

The latter form of the functional equation is related to the expression of (^) s (Q(i)(s) as the right 
Mellin transform of a tempered distribution which is supported in [0, oo) and which is self-reciprocal 
under the 7i-transform: 



(-) s (®(i)( s ) = / (T^ c m 6 nm (x) - -)x~"dx 

77 J ° rn>l 4 



(16a) 



r°° i i 

/ J (2y/xy)C^2 c «^«(y) ~ 4) d y = X! c ™^™( x ) - ^xxiix) = E(x) (x > 0) (16b) 

*' n>l m>l 

The invariance of E under the 7i-transform is equivalent to the validity of the functional equation 
of (^) s CQ(i)( s ) an d ^ having a pole with residue j at s = 1. The co-Poisson intertwining becomes 
the assertion: 

^ Jo(2V^y) I £ Cm — - 7 / 
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(17) 

If g is smooth with support in [b, B], < b < B < 00, then we have on the right hand side 
a function of Schwartz decrease at +00 (compare to Theorem 3), and its 7i-transform is also of 
Schwartz decrease at +00. The former is constant for < y < ttB^ 1 and the latter is constant for 
< x < irb. The supremum of the values obtainable for the product of the lengths of the intervals 
of constancy is 7r 2 . But, as for the cosine and sine transforms, there does exist smooth functions 
which are constant on a given (0, a) for arbitrary a > with an TL transform again constant on 
(0, a) and have Schwartz decrease at +00 (the two constants being arbitrarily prescribed.) 

De Branges and V. Rovnyak have obtained [1, 28] rather complete results in the study of the 
Hankel transform of order zero f(x) 1— > g{y) = / °° y/xyJo(xy)f(x)dx from the point of view of 
understanding the support property of being zero and with transform again zero in a given interval 
(0,6). They obtained an isometric expansion (Theorem 1 of section 2) and also the detailed de- 
scription of the related spaces of entire functions ([1]). The more complicated case of the Hankel 
transforms of non-zero integer orders was treated by J. and V. Rovnyak [29, 30]. These, rather 
complete, results are an indication that the Hankel transform of order zero or of integer order is 
easier to understand than the cosine or sine transforms, and that doing so thoroughly could be 
useful to better understand how to try to understand the cosine and sine transforms. 

The kernel Jo(2y / uI)) of the H-transform satisfies the Klein-Gordon equation in the variables 
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x = v — u, t = v + u: 

(j^T + !) J o(2V^) = (□ + l)Jo(2v^) = - ^ + l)J (Vi 2 -^ 2 ) = (18) 

It is a noteworthy fact that the support condition, initially considered by de Branges and V. Rovnyak, 
and which, nowadays, is also seen to be in relation with the co-Poisson identities, has turned out to 
be related to the relativistic causality governing the propagation of solutions to the Klein-Gordon 
equation. This has been established in [9] where we obtained as an application of this idea the 
isometric expansion of [1, 28] in a novel manner. It was furthermore proven in [9] that the H. 
transform is indeed an (absolute) scattering, in fact the scattering from the past boundary to the 
future boundary of the Rindler wedge < \t\ < x for solutions of a first order, two-component 
("Dirac"), form of the KG equation. 

In the present paper, which is completely independent from [9], we shall again study the H- 
transform and show in particular how to recover in yet a different way the earlier results of [1, 28] 
and also we shall extend them. This will be based on the methods from [5, 6], and uses the 
techniques motivated by the study of the co-Poisson idea [8] . Our exposition will thus give a fully 
detailed account of the material available in summarized form in [5, 6]. Then we proceed with a 
development of these methods to provide the elucidation of the (two dimensions bigger) spaces of 
functions constant in (0, a) and with 7i-transforms constant in (0, a). 

The use of tempered distributions is an important point of our approach 5 ; also one may envision 
the co-Poisson idea as asking not to completely identify a distribution with the linear functional 
it "is" . In this regard it is of note that the distributions which arise following the method of [5] 
are seen in the present case of the study of the 7Y-transform to have a very natural formulation as 
differences of boundary values of analytic functions, that is, as hyperfunctions [23]. We do not use 
the theory of hyperfunctions as such, but could not see how not to mention that this is what these 
distributions seem to be in a natural manner. 

The paper contains no number theory. And, the reader will need no prior knowledge of [2]; 
some familiarity with the to- function of Hermann Weyl [10, 21, 26] is necessary at one stage of 
the discussion (there is much common ground, in fact, between the properties of the to- function 
and the axioms of [2]). The reproducing kernel in any space with the axioms of [2] has a specific 
appearance (equation (109) below) which has been used as a guide to what we should be looking 
for. The validity of the formula is re-proven in the specific instance considered here 6 . Regarding 
the differential equations governing the deformation, with respect to the parameter a > 7 , of the 

5 at the bottom of page 456 of [1] the formulas given for A(a, z) and B(a, z) as completed Mellin transforms are 
lacking terms which would correspond to Dirac distributions; possibly related to this, the isometric expansion as 
presented in Theorem II of [1] is lacking corresponding terms. The exact isometric expansion appears in [28] and the 
exact formulas for A(a, z) and B(a, z) as completed Mellin transforms appear, in an equivalent form, in [30, cq. (37)]. 

e the critical line here plays the role of the real axis in [2], s is i — iz and the use of the variable s is most useful 
in distinguishing the right Mellin transforms which need to be completed by a Gamma factor from the left Mellin 
transforms of "theta"-like functions. 

7 the a here corresponds to ^a 2 in [1], 
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Hilbert spaces, we depart from the general formalism of [2] and obtain them in a canonical form, as 
defined in [21, §3]. Interestingly this is related to the fact that the A and B functions (connected 
to the reproducing kernel, equation (109)) which are obtained by the method of [5] turn out not 
to be normalized according to the rule in general use in [2]. Each rule of normalization has its 
own advantages; here the equations are obtained in the Schrodingcr and Dirac forms familiar from 
the spectral theory of linear second order differential equations [10, 21, 26]. This allows to make 
reference to the well-known Weyl-Stone-Titchmarsh-Kodaira theory [10, 21, 26], and to understand 
H as a scattering. Regarding spaces with the axioms of [2], the articles of Dym [14] and Rcmling 
[27] will be useful to the reader interested in second order linear differential equations. And we refer 
the reader with number theoretical interests to the recent papers of Lagarias [18, 19]. 

The author has been confronted with a dilemma: a substantial portion of the paper (most of 
chapters 5, 6, 8) has a general validity for operators having a kernel of the multiplicative type 
k(xy) possessing certain properties in common with the cosine, sine or H transforms. But on the 
other hand the (essentially) unique example where all quantities arising may be computed is the 
Ti. transform (and transforms derived from it, or closely related to it, as the Hankel transforms of 
integer orders). We have tried to give proofs whose generality is obvious, but we also made full use 
of distributions, as this allows to give to the quantities arising very natural expressions. Also we 
never hesitate using arguments of analyticity although for some topics (for example, some aspects 
involving certain integral equations and Fredholm determinants) this is certainly not really needed. 

2 Hardy spaces and the de B ranges- Rovnyak isometric ex- 
pansion 

Let us state the isometric expansion of [1, 28] regarding the square integrablc Hankel transforms of 
order zero. We reformulate the theorem to express it with the H transform (11) rather than the 
Hankel transform of order zero. 

Theorem 1 ([1], [28]). : Let k G L 2 (0, oo; dx). The functions fi and gi, defined as the following 
integrals: 




(19a) 



gi{y) = Kv) - \ J—^—Ji{ 2 Vy( x - v))K x ) dx , 



Jy v x ~y 



(19b) 



exist in L 2 



in the sense of mean-square convergence, and they verify: 




(19c) 



The function k is given in terms of the pair (/i, <?i) as: 



k(x)=g 1 (x) + J (2 \J y{x - y))fi{y) dy - 





x-y 



I! 



Ji(2v / y(x - y))gi(y) dy 



(19d) 
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The assignment k i— ► (/i, yi) is a unitary equivalence o/L 2 (0, oo; dx) £ 2 (0, oo; dy)©L 2 (0, oo; dy) 
smc/i that the TL-transform acts as (/i,<?i) i— > (fifi, Furthermore k and TL{k) both identically 
vanish in (0,a) ij and only if fi and gi both identically vanish in (0,a). 

Let us mention the following (which follows from the proof we have given of Thm. 1 in [9]): if 
fit fi> 9i> 9i are i n L 2 then k, k 1 and 7i(k)' are in L? . Conversely if k, k' and H(k)' are in L 2 then 
the integrals defining fx (y) and g\ (y) are convergent for each y > as improper Riemann integrals, 
and fi and g[ are in L 2 . 

It will prove convenient to work with (f(x),g(x)) = |(yi(§) + /i(f ),5i(§ ) - /i(f )): 

f(y) = + lf y °° 2 (joWv&c-v)) - ^^-MVy(2*-y))) Hx)dx (20a) 



(20b) 



9(V) = \k{\) - \ jH (joWvfr-v)) + ^^J—J^y/ypx-y)?) k(x)dx 

k{x) = f(2 X ) + \J^ (joiVypx-y)) - \j^—y J ^y^ x -y))) /(») d y 

+ 5 (2x)-i^ (joWvV* - y)) + \/^^ J ^y( 2x - y))) 9(y) d v ( 20c ) 

\k(x)\ 2 dx= / \f(y)\ 2 + \g(y)\ 2 dy (20d) 



The TL transform on k acts as (/, y) i— * (f,—g)- The pair (k,H(k)) identically vanishes on (0, a) 
if and only if the pair (/, g) identically vanishes on (0,2a). The structure of the formulas is more 
apparent after observing {x,y > 0): 

M\/y( 2x - y))io< y <2x(y)) = <My) - ^ y/ 2 s- y Jl ^ y ^ 2x ~ y)) 1 o< j y<2 3; (y) (21) 

In this section I shall prove the existence of an isometric expansion k «-* (/, g) having the stated 
support properties and relation to the 7i-transform; that this construction does give the equations 
(20a). (20b). (20c), will only be established in the last section (9) of the paper. The method 
followed in this section coincides partly with the one of V. Rovnyak [28] ; we try to produce the most 
direct arguments, using the commonly known facts on Hardy spaces. The reader only interested in 
Theorem 1 is invited after having read the present section to then jump directly to section 9 for the 
conclusion of the proof. 

To a function k £ L 2 (0, oo; dx) we associate the analytic function 

/>oo 

jfc(A) = / e lXx k(x)dx (3(A) > 0) (22) 



with boundary values for A 6 I again written A; (A), which defines an element of L 2 (R, |^), the 
assignment k i— > k being unitary from L 2 (0, oo; dx) onto H 2 (9(A > 0), 7p). Next we have the 
conformal equivalence and its associated unitary map from H 2 (9(A > 0), |£) to H 2 (|w| < 1, |^): 

w = ^ K{w) = ±±±±%(\) (23) 
A + 1 V2 i 



9 



It is well known that this indeed unitarily identifies the two Hardy spaces. With ko(x) = e 
fco(A) = jip-, K Q {w) = and ||fc || 2 = / °° er 2x dx = \ = ||A' || 2 - The functions fc„(A) = 
(xrt)"xTi corres P on d to K n {w) = -^w n . To obtain explicitely the orthogonal basis (fc n )n>0) we 
first observe that w = 1 — 2-^pj, so as a unitary operator it acts as: 

f-X C-X 

w ■ k(x) = k(x) -2 e- (x - y) k{y) dy = k(x) - e~ x 2 / e y k{y) dy (24) 
Jo Jo 

Writing k n (x) = P n {x)e~ x we thus obtain P n+ i(x) = P n {x) — 2 J Q P n {y) dy: 

So as is well-known P n (x) = Ln\2x) (in the notation of [31, §5]) where the Laguerre polynomials 
Ln \x) are an orthonormal system for the weight e~ x dx on (0, oo). 

One of the most common manner to be led to the 7i-transform is to define it from the two- 
dimensional Fourier transform as: 

H{f){\r 2 ) = ±-JJ e ^y^ x ^f( y -l±^)d yi dy 2 = (J*'*" "^) f(l s2 ) sds 

(26) 

n{f){-r 2 )=j^ J (rs)f(- S 2 )sd s r 2 = xl + x 2 2 ,s 2 = yj + y 2 2 

which proves its unitarity, self-adjointness, and self-reciprocal character and the fact that it has 
er x has self-reciprocal function. The direct verification of 7i(fco) = fco is immediate: Ti.(ko)(x) = 
/ °° J (2^y)e-y dy = £^° =0 ^fx n y» e -« dy = e~ x . Then, W(e" te ) = t^e"? for each t > 0. 
So J °° e~ tx H(k)(x) dx — t" 1 e~i x k(x) dx hence: 

Vfc G L 2 (0, oo; dx) H{k)(X) = ^-k(^-) (27) 

A A 

With the notation TC(K) for the function in H 2 (|u>| < 1) corresponding to TL(k), we obtain from 
(23), (27), an extremely simple result: 8 

H{K){w) = K{-w) (28) 
This obviously leads us to associate to K(w) — X)^Lo c n wTl the functions: 

oo 

F{w) := C2„w n (29a) 

n=0 

oo 

G(w) := C2n+iw n (29b) 

K(w) = F(w 2 ) + wG{w 2 ) (29c) 

and to k the functions / and g in £ 2 (0, +oo; dx) corresponding to F and G. Certainly, ||fc|| 2 = 
ll/ll 2 + llffl| 2 j an d the assignment of (f,g) to k is an isometric identification. Furthermore, certainly 



'we also take note of the operator identity H ■ w = — w ■ Ti. 
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the Tl transform acts in this picture as (/, g) t— ► (/, —g). Let us now check the support properties. 
Let a(m) be the leftmost point of the (essential) support of a given m <E L 2 (0,oo;dx). As is 
well-known, 

— a(m) = limsup — log |m(i£)| , (30) 

t— >+oo t 

If w corresponds to A via (23) then w 2 corresponds to ^(A — y), so if to a function / with corre- 
sponding F(w) we associate the function tp(f) G L 2 (0, oo; dx) which corresponds to F(w 2 ), 

(t + t) = (^±i + i±i) , (si) 

then we have the identity: 

= (32) 

Returning to F (rcsp. /) and G (resp. g) associated via (29a). (29b), to K (resp. k) we thus have 
k = i-'if) + w ' VKff); W(fc) = ^(/) — W ■ i[>(g), hence if the pair (/,<?) vanishes on (0,2a) then the 
pair (k,7i(k)) vanishes on (0, a) (clearly the unitary operator of multiplication by w = yqj| does not 
affect a(m).) Conversely, as a(f) = 2a(k + H(k)) and a(g) = 2a(k — H(k)), if the pair (k, TL(k)) 
vanishes on (0,a) then the pair (/, g) vanishes on (0,2a). 

We have thus established the existence of an isometric expansion, its support properties, and its 
relation to the 7i-transform. That there is indeed compatibility of (20a) and (20b) with (29a) and 
(29b), and with (20c), will be established in the last section (9) of the paper with a direct study of 
(31). In the meantime equations (20a), (20b), (20c) and (20d) will have been confirmed in another 
manner. Yet another proof of the isometric expansion has been given in [9]. 



3 Tempered distributions and their TL and Mellin transforms 

Any distribution D on R has a primitive. If the closed support of D is included in [0, +oo), then 
it has a unique primitive, which wc will denote D(x) dx, or, more safely, Z?( -1 ), which also has 
its support in [0, +oo). The temperedness of such a D is equivalent to the fact that D(~ N > for 
N ^> is a continuous function with polynomial growth. With D^-~ N ^(x) = (1 + x 2 ) M g(N,M)( x )> 
M»0, we can express D as P(x, -£)(g) where P is a polynomial and g £ L 2 (0, oo; dx). Conversely 
any such expression is a tempered distribution vanishing in (— oo,0). The Fourier transforms of 
such tempered distributions D(X) appear thus as the boundary values of Q(^,A)/(A) where Q 
are polynomials and the f's belong to H 2 (3(A) > 0). As taking primitives is allowed we know 
without further ado that this class of analytic functions is the same thing as the space of functions 
g(X) = R(j%, A, A _1 )/(A), R a polynomial and / € H 2 . It is thus clearly left stable by the operation: 

g »H(g)(A):=jg(=±) (3(A) > 0) (33) 
which will serve to define the action of Ti. on tempered distributions with support in [0, +oo). 
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Let us also use (33), where now A S i, to define H as a unitary operator on L 2 (— oo, +00; dx). 
It will anti-commute with f(x) — > /(— x) so: 

/■OO 

«(/)(*) = / (Jo(2v^)lx>o(!c)lv>o(l/)-Jb(2V5i0lx<o(x)l tf <o(»))/(tf)di/ (34) 



Useful operator identities are easily established from (33): 

x— ■ H = —Ti. ■ —x and — x ■ TL = —Ti. ■ x— (35a) 

ax ax dx dx 

d f x f x d 

— ■H = H ■ / and -H = H ■ — (35b) 

dx ,/n Jo dx 



x ■ Ti. = —Ti ■ — x — and TC ■ x — — — x — ■ Ti (35c) 

dx dx dx dx 

It is important that ^ is always taken in the distribution sense. It would actually be possible 

to define the action of Ti on distributions supported in [0, +oo) without mention of the Fourier 

transform, because these identities uniquely determine Ti(D) if D is written {-£:) N (1 + x) M gN ,m(x) 

with gN,M £ L 2 (0, oo; dx). But the proof needs some organizing then as it is necessary to check 

independence from the choice of TV and M, and also to establish afterwards all identities above. So 

(33) provides the easiest road. Still, in this context, let us mention the following which relates to 

the restriction of H(D) to (0, +oo): 

Lemma 2. Let k be smooth on R with compact support in [0, +oo). Then TL{k) is the restriction 
to [0, +oo) of an entire function 7 which has Schwartz decrease as x — > +00. For any tempered 
distribution D with support in [0, +00), there holds 



H{D)(x)k{x) dx = / D(x)~f{x) dx , (36) 
Jo 

where in the right hand side in fact one has D(x)j(x)Q(x) dx where the smooth function is 1 
for x > — I and for x < — | and is otherwise arbitrary (as is e). 

Let us suppose k = for x > B. Defining: 

7 (x)= / Jo(2^)k(y)dy (37) 
Jo 

we obtain an entire function and, according to our definitions, TL{k){x) = ~f(x)\ x> o{x) as a distribu- 
tion or a square-integrable function. Using (35c) {TL = —\Ti- 4- x ^ for x > 0) and bounding J by 
1 we sec (induction) that 7 is 0(x~ N ) for any N as x — > +00 , and using (35a) (-fr • Ti. = —-^H • -j^x 
for x > 0) the same applies to its derivative and also to its higher derivatives. So it is of the Schwartz 
class for x — > +00. 

Replacing D by TC(D) in (36) it will be more convenient to prove: 

/•OO 

D{x)k{x) dx = I H(D){x)~f{x)dx (38) 
Jo 
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If (38) holds for D (and all fc's) then < D',k>= - < D, k' > = - < H(D), -6{x) $™ j(y)dy > (ob- 
serve that f* j(y) dy = H(k')(x) vanishes at +00) so < D',k >= + < f*H(D),9~f >=<H(D'),9~f > 
hence (38) holds as well for D' (and all k's). So we may assume D to be a continuous function of 
polynomial growth. It is also checked using (35c) that if (38) holds for D it holds for xD. So we 
may reduce to D being square-intcgrable. and the statement then follows from the self-adjointncss 
of TL on L 2 (or we reduce to Fubini). 

The behavior of TL with respect to the translations r a : f(x) 1— > f(x — a) is important. For 
/ G L 2 (M; dx) the value of a is arbitrary and we can define 

t* :=TLr a TL (39a) 
7^7)(A)=e <oA /(A) (39b) 
T*(f)(\)=e™^f(\) (39c) 

We observe the remarkable commutation relations (which would fail for the cosine or sine trans- 
forms): 

Va,6 Ta T* = T*T a (40) 

For a distribution D the action of rf is here defined only for a > —a(TL(D)), where a(E) is the 
leftmost point of the closed support of the distribution E. On this topic from the validity of (30) 
when / £ L 2 (0, 00; dx), and invariance of a under derivation 9 , integration, and multiplication by x, 
one has: 

-a(E) = limsup - \og\E(it)\ (41) 

i->+oo t 

We thus have the property, not shared by the cosine or sine transforms: 

a > —a(H(D)) a(r#(D)) = a(D) (42) 

We now consider D with a(D) > and a(TC(D)) > and prove that its Mcllin transform is 
an entire function with trivial zeros at 0, —1, —2, following the method of regularization by 
multiplicative convolution and co-Poisson intertwining from [8]. The other, very classical in spirit, 
proof shall be presented later. The latter method is shorter but the former provides complementary 
information. 

In [8, §4. A] the detailed explanations relative to the notion of multiplicative convolution are 
given: 

(9*D)(x) " = " feWDfy* (43) 
where we will in fact always take g to have compact support in (0, +00). It is observed that 

g*xD = x(-*D) (g*D)' = -*D' (44) 

x x 



9 It is important in order to avoid a possible confusion to insist on the fact that ^ is always taken in the distribution 
sense so for example -^~S- X >0 = $( x ) indeed has the leftmost point of its support not affected by -4-. 
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The notion of right Mellin transform J °° D(x)x s dx is developed in [8, §4.C], for D with support 
in [a, +00), a > 0: 

D(s) = s(s + 1) • • • (s + JV - l)5<r-W(s + N) , (45) 

where N 3> 0. The meaning of D is as the maximal possible analytic continuation to a half-plane 
3?(s) > a, where a is as to the left as is possible. The notion is extended 10 in [8, §4.F] to the case 
where the restriction of D to (—a, a) is "quasi-homogeneous". For example, if D\,_ a a ^ — lo<a:<a 
(resp. S), then D is defined as D\ with D\ = D — lo<a:<oo (resp. D — 8.) Then, also in the extended 
case, the following holds: 

^Td(s) = g(s)D(s) (46) 

where g in an integrable function with compact support in (0, 00) and g{s) is the entire function 
Jo°° 9(t)t s We then have the following theorem: 

Theorem 3. Let D a tempered distribution with support in [a, +00), a > and such that Ti(D) 
also has a positive leftmost point of support. Let g be a smooth function with compact support in 
(0,oo). Then the multiplicative convolution g * D belongs to the Schwartz class. 

This is the analog of [8, Thm 4.29]. The function k(t) = (Ig)(t) = i s defined and it is 

written as k = H(7l a:> o) where 7 is the entire function, of Schwartz decrease at +00 such that 
TC(k) = 7 • l x >o- Then it is observed that 

t>0 (g*D)(t)= [ D{x)^^dx= ( H(D)(x)j(tx)dx (47) 

Jo t Jo 

We have used Lemma 2. Then the Schwartz decrease of J °° 7i(D)(x)^f(tx) dx as t — > +00 is es- 
tablished as is done at the end of the proof of [8, Thm 4.29], integrating by parts enough times to 
transform Tt(D) into a continuous function of polynomial growth, identically zero on [0, c], c > 0. 

Theorem 4. Let D a tempered distribution with a positive lefmost point of support and such that 
Ti(D) also has a positive leftmost point of support. Then D(s) and T(s)D(s) are entire functions 
and: 

T(s)D(s)=T(l-s)n{D)(l-s) (48) 

We first establish: 

Theorem 5 ("co-Poisson intertwining"). Let D be a tempered distribution supported in [0, +00) 



and let g be an integrable function with compact support in (0,oo). Then, with (Ig)(t) 



1 



H(g * D) = {Lg) * H(D) (49) 



10 if D is near the origin a function with an analytic character, then straightforward elementary arguments allow a 
complementary discussion. However if D is just an element of L 2 (0, 00; dx) then D is a square-integrable function on 
the critical line, and nothing more nor less. 
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Let us first suppose that D is an L 2 function. In that case, we will use the Mellin-Plancherel 
transform / i— ► f(s) = J °° f(t)t~ s dt, for / square integrable and 5R(s) = |. Then g * f is, changing 
variables, the Fourier transform of an additive convolution where one of the two has compact support, 
well known to be the product g ■ f. We need also to understand the Mellin transform of H(f). Let 
us suppose ft(x) = cxp(-tx). Then TL(fi) = 7/1 has Mellin transform H(ft)(s) — t~ s T(l — s) and 
ft(s) = i s_1 r(l — s), so we have the identity for such f's: 

nf)(s) = ^fAf(i-s) (50) 

r(s) 

The linear combinations of the ft's are dense in L 2 , so (49) holds for all f's as an identity of 
square integrable functions on the critical line. We are now in a position to check the intertwining: 
K(g~7})(s) = ^fg(l - s)f(l - s) = hj(a)H(JM = Ig^H{f){s). 

For the case of an arbitrary distribution it will then be sufficient to check that if (49) holds for 
D it holds for xD and for D' . This is easily done using (44). We have g * (£)') = (xg * D)' , so 
H{g*D')= f*H(xg*D) = f*(^*H(D)) = Ig*(f* H(D)) = Ig*H(D'). A similar proof is done 
for xD. This completes the proof of the intertwining. 

The theorem 4 is then established as is [8, Thm 4.30]. We pick an arbitrary g smooth with 
compact support in (0, 00). We know by theorem 3 that g * D is a Schwartz function as x — > +00, 
and certainly it vanishes identically in a neighborhood of the origin, so g * D(s) = g(s)D(s) is an 
entire function. So D(s) is a meromorphic function in the entire complex plane, in fact an entire 
function as g is arbitrary. We then use the intertwining and (50) for square integrable functions. 
This gives g(l - s)H(D)(s) = Ig*H(D)(s) = H{g*~D)(s) = ^^-g{\ - s)D(l - s). Hence, indeed, 
after replacing s by 1 — s: 

T(s)D(s)=T(l-s)H(D)(l-s) (51) 

The left-hand side may have poles only at 0, —1, . . . , and the right-hand side only at 1, 2, ... . So 
both sides are entire functions and D(s) has trivial zeros at 0, —1, —2, . . . 

We now give another proof of Theorem 4, which is more classical, as it is the descendant of the 
second of Riemann's proof, and is the familiar one from the theory of theory of L-serics and modular 
functions. The existence of two complementary proofs is instructive, as it helps to better understand 
the role of the right Mellin transform J °° f(x)x~ s dx vs. the left Mellin transform #(«i)i s_1 dt. 

To the distribution D we associate its "theta" function 11 #d(A) = -D(A) = / °° e lXx D(x) dx, 
which is an analytic function for SJ(A) > 12 . Right from the beginning we have: 

0H(D)(it) = (52) 

If the leftmost point of the support of D is positive then 6d {it) has exponential decrease as t — > +00 
and 0£)(rf)t s_1 dt is an entire function. If also the leftmost point of support of TL{D) is positive 

11 the author hopes to be forgiven this temporary terminology in a situation where only the behavior under A 1— » 
is at work. 

12 we adopt the usual notation, and consider djj as a function of it rather than t. 
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then 0-H(£>)(it) nas exponential decrease as t — > +00 and Jq 1 #D(ii)i s_1 di = /j 30 9i-i(D){it)t~ 3 dt is an 
entire function. So, under the support property considered in Theorem 4 T>(s) := J °° $D(it)t s_1 dt 
is indeed an entire function, and the functional equation is 

V{s) = V*(l - s) (53) 

with v*( s ) = / °° e H{D) {t)t s - 1 dt. 

To conclude we also need to establish: 

V{s) = T(s)D(s) (54) 

We shall prove this for 5R(s) 3> under the hypothesis that D has support in [a, +00), a > (no 
hypothesis on H(D)). In that case, as 9d(H) is 0(tr N ) for a certain TV as i — * (t > 0), and is of 
exponential decrease as t — > +00, we can define 2?(s) = / °° #£)(it)i s_1 df as an analytic function for 
3?(s) 3> 0. Let us suppose that D is a continuous function which is 0(x~ 2 ) as x — > +00. Then, for, 
3?(s) > 0, the identity (54) holds as an application of the Fubini theorem. We then apply our usual 
method to check that if (54) holds for D it also holds for xD and for D' . For this, obviously we 
need things such as D'(s) = sD(s + 1) [8, 4.15] and xD(s) = D(s — 1), the formulas 6d' = — i^@D, 
9 X £i = —i-7^9]j, and T(s + 1) = sr(s). The verifications are then straightforward. 

In summary we have seen how the support property for D and Ti-(D) is related in two comple- 
mentary manners to the functional equation, one using the right Mellin transform D(s) of D and 
the idea of co-Poisson, the other using the left Mellin transform T>(s) of the "theta" function Op 
associated to D as an analytic function on the upper half- plane and the behavior of 0jj{it) under 
i h> 7. It is possible to push further the analysis and to characterize the class of entire functions 
T>(s) = T(s)D(s), as has been done in [8] in the case of the cosine and sine transforms. It is also 
explained in [8] how the discussion extends to allow finitely many poles. The proofs and statements 
given there are easily adapted to the case of the Tt transform. Only the case of poles at 1 and 
will be needed here and this corresponds, either to the condition that D and Ti-(D) both restrict in 
(—a, a) for some a > to multiples of the Dirac delta function, or, that they are both constant in 
[0, a) for some a > 0. We recall that the Mellin transform D(s) is defined in such a manner, that it 
is not affected from either substracting S or l x >o from D. 

4 A group of distributions and related integral formulas 

We now derive some integral identities which will prove central. The identities will be re-obtained 
later as the outcome of a less direct path. We are interested in the tempered distribution g a (x) 
whose Fourier transform is exp(ia^). Indeed rf{f) (equation (39a)) is the additive convolution of 
/ with g a : we note that g a differs from S(x) by a square intcgrablc function as 1 — cxp(— iaX^ 1 ) = 
0|A|-foo(|^| _1 ); so there is a convolution formula rf(j) = f — f a * f for a certain square integrablc 
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function f a . For / G L 2 , the convolution f a * f as the Fourier transform of an Z^-function is 
continuous on R. Starting from the identity exp(ia^-) = — i\j exp(ia^-) we identify g a for a > 
as -^TL8 a . It is important that ^ is taken in the distribution sense. So we have, simply: 



a J\ (2^/ax) 
y/ax 

If a < then g a (x) = g- a (-x), f a (x) = f- a (-x). So 



g a (x)=5(x)- ulK -^"' l x>0 (x) (a>0) (55) 
\f ax 



g. a (x) = 5(x) - aJl{2 / f^ ) l x<0 {x) {-a < 0) (56) 

V — ax 

The group property under the additive convolution g a * g = g a +b leads to remarkable integral 
identities f a +b = fa + fb — fa* fb involving the Bessel functions. The pointwise validity is guaranteed 
by continuity; the Planchcrcl identity confirms the identity, where f a {%) — ~ Jl ^j a ^^ ~ ^-x> (%) for 
a > and /_ a (:r) = f a (-x): 

fa+b = fa + fb - fa * fb (57) 

At x = the pointwise identity is obtained by continuity from either x > or x < 0. We have 
essentially two cases: g a * gb for a, 6 > and g a * o_b for a > b > 0. The following is obtained: 

Proposition 6. Let a > b > and x > 0. There holds: 



{a + b).h{2^/(a + b)x) _ a J^y/Ex) , &Ji(2\/&c) /"* aJx(2^/ay) bJ^y/Hx - y)\ ^ 



y/(a + b)x Vox y/ax Jo V®y V b (x - 2/) 



(a- 6)^(2^/(0 -b)a;) _ aJtfy/ZZ) f°° aJ^^ay) bJ 1 (2 x /b(y - x)) 



y/(a-b)x Vox J x y/ay ^fb{y - x) 



Q = bJtf-Jbx) _ f 00 aJi(2^) 6Jx(V^+g) 

Exchanging a and b and changing variables we combine (58b) and (58c) into one single equation 
for x > and a, 6 > 0: 



(a-b)J 1 (2 v /(a-b) a; ) i aJ 1 (2^i) aJi(2y/a(y + g) frJi(2y^) 

l a _ b >o(a-6) = = / = — dy (59) 

V (a - b)x V ax Jo \/a{y + x) Vby 



The formula for x = in (59) is obtained by continuity. It is equivalent to 

d 1 1 

Ji(u)Ji(cu) — = - min(c, -) (c> 0) (60) 
u 2 c 

which is a very special case of formulas of Weber, Sonine and Schafhcitlin ([33, 13.42.(1)]). Another 
interesting special case of (59) is for a = b. The formula becomes 

Ji(2Vz) _ f°° Ji(2V») Ji(2V^+y) ^ 



which is equivalent to a special case of a formula of Sonine ([33, 13.48.(12)]) 
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We already mentioned the equation Jo(2y/uv) = —Jq(2^/uv). New identities are obtained 
from (59) or (58a) after taking cither the a or the b derivative. We investigate no further (59) as the 
corresponding semi-convergent integrals, in a form or another, are certainly among the formulas of 
[33, §13]. Let us rather focus more closely on the case a, b > ((58a).) We have a function which 
is entire in a, b, and x and the identity holds for all complex values of a, b, and x. Let us take the 
derivative with respect to a: 



J (2V(a + b)x) = J (2V^) - / JoVVw) bJl ^ b( - X y) l dy (62) 

Jo x /b{x-y) 



We replace b by —b and then set x = b. This gives: 



J (2 VKb^)) = J (2Vba~) + [ J (2^Ey) ^^—J^ dy (63) 

Jo VHb-y) 



We take the derivative of (62) with respect to b: 



= - / J Q {2^Elj)J {2y / b{x - y)) dy 

^/{a + b)x Jo 



Then we replace b by —b and set x = b: 



y'b(b-a) Jo 



(64) 



(65) 



Combining (63) and (65) by addition and substraction we discover that we have solved certain 
integral equations: 



= Io(2VHb=xj) bh{ ^t_ x)) = (1 + A) Jo ( 2v /&(&^)) (66a) 



^(x) = / (2v^(6^)) + bIli ^— X)) = (1 - A )/o(2v ^(^)) (66b) 



#(z) + / Jo{2^y)4>+(y) dy = J (2Vbx) (66c) 



n 



<f>b( x )~ M2y/xy)^{y)dy = J (2Vbx) (66d) 
Jo 

The significance will appear later in the paper and we leave the matter here. The method was devised 
after the importance of solving equations (66c) and (66d) had emerged and after the solutions (66a) 
and (66b) had been obtained as the outcome of a more indirect path. Of course, direct verification 
by replacement of the Besscl functions by their series expansions is possible and easy. 

5 Orthogonal projections and Hilbert space evaluators 

Let a > and let P a be the orthogonal projection on L 2 (0, a; dx) and Q a = TCP a TC the orthogonal 
projection on H.(L 2 (0,a;dx)) and let K a C L 2 (0, oo; dx) be the Hilbert space of square integrablc 
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functions / such that both / and H{f) have their supports in [a, oo). Also we shall write H a = 
P a HP a . Also we shall very often use D a = Hi = P a HP a HP a . Using: 

Jo(2^y) = f>l)"^f , (67) 

we exhibit H a = PaTi.P a as a limit in operator norm of finite rank operators so P a Ti.P a is a compact 
(self-adjoint) operator. It is not possible for a non zero / G L 2 (0,a;dx) to be such that ||if (/)|| = 
11/11, as this would imply that H a (f) vanishes identically for x > a, but H a (f) is an entire function. 
So the operator norm of H a is strictly less than one, and 1 ± H a as well as 1 — D a arc invcrtiblc. 
We consider the equation 

<j> = u + H{v) u,v G L 2 (0,a;dx) (68) 

Hence: 

u + H a (v)=P a (4>) (69a) 
H a (u)+v = P a {H{4>)) (69b) 
u={l-D a )~ 1 (P a {^)-H a P a H{4>)) (69c) 

V = (1 - Da^i-HaPa^) + P a H{^)) (69d) 

Then if <f> n = u n + H(v n ) is L 2 -convergent, (u n ) and (v n ) will be convergent, and the vector space 
sum L 2 (0, a; dx) +H(L 2 (0, a; dx)) is closed. Its elements are analytic functions for x > a so certainly 
this is a proper subspace of L 2 . Hence we obtain that each K a is not reduced to {0} and 

= L 2 (0, a; dx) + TL(L 2 {Q, a; dx)) (70) 

We also mention that \J a >oK a is dense in but not equal to L 2 (0, oo; dx), more generally that U a> bK a 
is dense in but not equal to K^, and also obviously n a<00 A' a = {0}, n a< bK a = Kj,. 

In this section a > will be fixed (all defined quantities and functions will depend on a, but 
this will not always be explicitely indicated.) We shall be occupied with understanding the vectors 
I g a G K a such that 



V/GK a / f(x)X:(x)dx = f(s)= f(x)x- s dx (71) 

J a J a 

and in particular we are interested in computing 

X a (s,z)= X a s {x)X*{x)dx (72) 



As a is fixed here, we shall drop the superscript a to lighten the notation. For the time being we 
shall restrict to SR(s) > \ and we define X s to be the orthogonal projection to K a of l x>a (x)x~ s . 
As a preliminary to this study we need to say a few words regarding: 

/•OO 

g s (x) :=H{l x>a (x)x- s )= J {2^)y- S dy (73) 
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The integral is absolutely convergent for 3?(s) > §, semi- convergent for 3?(s) > ^, and g s is defined 
by the equation as an L 2 function for 0?(s) > i (it will prove to be entire in s for each x > 0). We 
need the following identity, which shows also that g s (x) is analytic in x > 0: 

/■a 00 n n+1 — s 

^Or)^)* 8 - 1 - / J^y/Zfiy-* dy = X (s)x s -' (74) 

This is obtained first in the range | < 5R(s) < 1: f°° Jo(2y/xy)y~ s dy = x s_1 J™ Jo(2^/y)y~ s dy = 
X s ' 1 (x(s) - Jo" M2VV)y~ s dy) = ^ _1 X(«) - £ Jo(2^xy)y- s dy. The poles at = 1, s = 2, 
. . . are only apparent. The identity is valid by analytic continuation in the entire plane 3?(s) > \. 
For each given x > we have in fact an entire function of s G C. But we are here more interested 
in g s as a function of x and we indeed see that it is analytic in C\] — 00, 0] (it is an entire function 
of x if s e — N). 13 

There are unique vectors u s , v s in L 2 (0, a; cix) such that 

l 2; >a(a;)a;- s = X s {x) + u s {x) + H(v s )(x) (75) 

and they are the solutions to the system of equations: 

u s + H a {v s ) = (76a) 
H a (u s ) + v s = P a (g s ) (76b) 

From (76a) we see that u s is in fact the restriction to (0, a) of an entire funtion, and from (76b) 
that v s is the restriction to (0,a) of a function which is analytic in C\] — 00, 0]. Redefining u s and 
v s to now refer to these analytic functions their defining equations become (on (0, +00)): 

u s +HP a (v s ) = (77a) 
HP a (u s ) + v s =g 8 (77b) 

and (75) becomes (we set X s (a) = X s (a+j): 

l x > a (x)x- s = X s (x) + l 0<x<a (x)u s (x) + HP a {v s ){x) (78a) 
l x > a {x)x~ 8 = X s (x) - l x > a (x)u s (x) (78b) 
X s (x) = l x > a (x)(x- s + u s (x)) (78c) 

The key to the next steps will be the idea to investigate the distribution (x-^ + s)X s on the (positive) 
real line. Let D s be x-^ + s. There holds: 

DM = -7Wi- B (79) 

To compute -^Pa{v s ) we first suppose 0?(s) > 1, so (we know the behavior as x — > from (74)) 
j^P a {v s ) = P a (v' s ) - v s (a)S a (x) and x-^P a (v s ) = P a (xv' s ) - av s (a)S a (x). This remains true for 



13 For some other transforms k(xy), such as the cosine transform, the argument must be slightly modified in order 
to accomodate the fact J^ 00 k(y)y~ s dy has no range of absolute convergence. 
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3?(s) > \. Applying D s to (77a) thus gives D s {u s ) — H(P a Di_ s (v s ) — av s (a)S a (x)). We similarly 
apply -Di-s to (77b) and obtain the following system: 



D s (u s )(x) - (nP a D 1 _ s v s )(x) = -av s (a)J (2y/ax) (80a) 
-(HP a D s u s )(x) + D^ s {v s ){x) = (D^ s g s ){x) - au s {a)J {2^i) (80b) 

From (73), we have D\- S g s = — HD s (l x>a x~ s ) = —H.(a 1 ~ s 5 a (x)) = — a 1 ^ 3 Jo(2y/ax). Let us define 

J*(x) = J (2y^x~) (81) 

We have proven: 

+D 3 u s - HP a Dx- s v s = -av s (a)J$ (82a) 
-HP a D s u s + D^ s v s = -a(a- s + u s (a))Jg (82b) 

Restricting to the interval (0, a) and solving, we find: 

P a D s u s = -o(l - £> Q )" 1 (w s (a) J a + (a~ s + u s {a))H a J*) (83a) 
P a D x _ s v s = -o(l - Da)- 1 ^- 8 + u s {a)) J a + v s (a)H a Jg) (83b) 

It is advantageous at this stage to define <j>£ and </>~ to be the solutions of the equations (in 
L 2 {0,a;dxj): 

4>t + H a 4>t = J% (84a) 

0a - Hata = ^0 (84b) 

We already know from (66a) and (66b) exactly what 4>t and 4>a are (in this special case of the Tt 
transform), but we shall proceed as if we didn't. We see from (84a), (84b) that 0+ and (j>~ arc 
entire functions, and we can rewrite the system as: 14 

<j>t +HP a 4>t = J% (85a) 
0- - HPacp- = J a (85b) 



We observe the identities: 



-1 ja 



{l-D a )- l J«=P a ^ a 2 ya (86a) 
(1 - D a )- l H a J- = P a 4,(1 (86b) 



So (83a) and (83b) become 



a s + u s (a) - v s (a) + a s + u s (a) + v s (a) 
D s u s = +a 4>Z - a 4> a (87a) 

aT s +u s (a) -v s (a) + q~ a + u s (a) + ^ s (q) 
£>i_ s w s = -a 4>J - a <p a (87b) 



14 in conformity with our conventions, these are identities on (0, oo); to see them as identities on C one must read 
J a Jo{2 s /xy)(j>t (y) dy rather than {HP a <t>t){x). 
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From (87a) we compute successively (again, these are identities on (0,+oo)): 

n-iP n a ~ s + u s{a) ~ Vs{a) , ia aT s + u s {a) + v s {a) 
HP a D s u s = a (J - 4>Z) - a (-J + a ) 

P a D s u s = a (d a -n(f>a) - a (-6 +rl<p a ) (89) 

In (89), Tt4>a should perhaps be more precisely written as Ti(4>^ l x yo). From (85a) we know that 
4>^l x >o is tempered as a distribution. From (78c) we compute D S X S = l x>a D s (u s ) + a(a~ s + 
u s (a))5 a (x) = D s u s — P a D s u s + a(a~ 3 + u s (a))8 a (x). From (87a) and(89) then follows: 

n y , a~ s + u s {a) - v s (a) + + a~ s + u s (a) + v s (a) 

D S X S = +a {4>l + H4>1 - S a ) - a - (<j> a - Ti(j> a + 6 ) , 

+a(a~ s + u s (a))S a (x) 

And the result of the computation is: 

n y , a~ s + u s ( a )-Vs(a) , A+ . . a~ s + u s {a) + v s (a) 

D S X S = +a (</>+ + H(j)Z) + a (-<f> a + H<j> a ) (91) 

We then define the remarkable distributions: 

4, = ^(#+tt#) (92a) 

-iB a = ^{-4,-+U4>-) (92b) 
E a = A a - iB a (92c) 

From (84a) we observe that A a has its support in [a, oo). Furthermore it is H. invariant. Similarly, 
—iB a , which is 7i anti invariant, also has its support in [a, +oo). We recover A a and —iB a from E a 
through taking the invariant and anti- invariant parts. We may also rewrite D S X S as: 

D S X S = V^(a- S + u s {a))E a - Vav s {a)HE a (93) 

Some other manners of writing A a and —iB a are useful: from (85a) Ti.^ = 5 a — P a 4>t an d from 
(85b) Ufa =5 a +P a 4>-, so: 

A a = ^(S a + <t>+l x>a ) (94a) 
-iB a = ^{5 a -ct>-l x>a ) (94b) 
And, we take also notice of the following definitions and identities: 

ja = Va(S a - 4>t 1 o<x<a) ja = ^faH<t>+ A a = i (j a + Hj a ) (95a) 

-ik a = y/a(S a + <f>a~\-Q<x<a) —ik a = \faH(j)~ B a = i (k a - Hk a ) (95b) 
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From (85a) and (85b) we know that </>+ and 4> a are bounded, so the right Mellin transforms are 
defined directly for 3i(s) > 1 by: 15 

Ms) = ^ (a- s + £ ri(x)x- s dx^j (96a) 

-W a (s) = ^ (a- s - £ 0~(x)x- s dx) (96b) 

Ea(s) = (a" s + i J y>+(x) ~ r a {x))x- s dx) (96c) 



_ i r°° 

H(E a )(s) = V^2 J (<t>i(x)+4>-(x))x- s dx (96d) 

From (85a) and (85b) we know that 4>^ — 0~ is square- intcgrable at +oo. so, using H.(p^ — 5 a — P a 4>i 
and Ti-4>a = + Pa^a wc compute: 



(4>t{x)-fc{x))x- s dx= I (H<t>i-Hct>-)g s (x)dx = - I (<£+(x) + 0" (*))$. (s) dx (97) 
Then using (86a): 



a ^M±^ (a:) g s (x) = - Z^)"Vo = f J Q (x)((l - D a )~ 1 (g s ))(x) dx 



JO 

(98) 

Comparing with (76a) and (76b) the right-most term of (98) may be written as J" Jo(2y/ax)v s (x) dx 
which in turn we recognize from (77a) to be — u s (a). We have thus proven the identity: 

Ea(s) = Va~(a- s +u s (a)) (99) 

In a similar manner we have: 

00 # (*> + *-<*w= r M 2v^)x-dx + r-rtto + +*to g .( X )*c (ioo) 



2 Ja 

a -<t>i{x)+j)-{x) 



(x)dx = I ((1-D a ) l H a J%)(x)g s (x)dx = - J J a (2y/ax)u s (x) dx 
= v s (a) — g s (a) — v s (a) — / Jo(2y /r ax)x~ s dx 

J a 

(101) 

A a (s) + i'Ba(s) = H(E a )(s) = yfc H t*M±4*W x -' dx = ^ Vs{a) (102 ) 



Then, we obtain the reformulation of (93) as: 

D S X S = E a (s)E a - H{EZ){s)HE a (103) 



15 the integral for E a (s) is certainly absolutely convergent for 5R(s) > ^ as <j>a — 4> a is square integrable on (0, oo), 
and in fact it is absolutely convergent for 5R(s) > j. As we know already completely explicitely (j>a and <f>a , we do not 
pause on this here. A general argument suitable to establish in more general cases absolute convergence for > cr 
for some a < = will be given later. 
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And, noting D s X s (z) = (s + z — l)X s (z) = (s + z — 1) X s (x)X z (x) dx we are finally led to the 
remarkable result: 

*„(.,*) = r x^x^dx = ^m-mnsmiz) m 

J a S + Z-1 

This equation has been proven under the assumption 5R(s) > 1, and ^St(z) > To complete the 
discussion we need to know that the evaluators / h- > /(s), s G C arc indeed continuous linear forms on 
K a . For 5R(s) > \, we have /(s) = J a °° /(x)ar s dx. For 5R(s) < \ we have /(s) = ^fc^-7Tf(7)(l — s). 
For 0?(s) = | continuity follows by the Banach-Steinhaus theorem, and of course more elementary 
proofs exist (as in [3] for the cosine or sine transform). So we do have unique Hilbcrt space vectors 
X a s G K a such that V/ G K a \fs G C f(s) = J a °° X^(x)f(x) dx. Then (104) holds throughout C x C 
by analytic continuation. 

The vectors X£ are zero for s G — N, and it is more precise to use vectors X" = T(s)X^ which 
are non-zero for all s G C. These vectors arc the evaluators 16 for / i— > J~(s), J-(s) = T(s)f(s). 
We recapitulate some of the results in the following theorem, whose analog for the cosine (or sine) 
transform was given in [5] (up to changes of variables and notations, the first paragraph as well as 
equation (108) are theorems from [1]; the equations (105), (106), (107) are our contributions. In 
this specific case of Ti we shall later identify exactly 4>~t an( i 0a anc ^ E a and £ a . As we shall explain 
the analog of the £ a -function in [1] has value 1 at s = i, and is not identical with the £ a here): 

Theorem 7. For a given a > let K a be the Hilbert space of square integrable functions f{x) on 
[a, +oo) whose Tt -transforms Jo(2y/xy)f(y)dy (in the I? '-sense) again vanish for < x < a. 
The completed right Mellin transforms T{s) f{s) = T(s) f{x)x~ s dx are entire functions and 
evaluations at s G C are continuous linear forms. 

Let Xg for each s G C be the unique vector in K a such that V/ G K a T(s) f{s) = J a °° f(x)X^(x) dx. 
Let and <p~ be the entire functions which are the solutions to: 



<t>tix) + / J o (2V^)0I(y) dy = J (2^) (105) 
Jo 

4>-(x)- [ M2^y-)4>-(y)dy = M2^) (106) 



o 



Then 



E a {s) = V^(a- S + \ ptetV) ' K{x))x- S dx) (107) 
is an entire function with trivial zeros at — N and, defining £ a {s) = T(s)E a (s), we have: 

Vfi, Z G C ^ X°(X)X°(X) dx = £ a{s)£a{z)-£a{l-s)£ a {l-z) ^ 
J a S -{- Z 1 

We knew in advance that we had to end up with a formula such as (108) (with a £ function to 
be discovered 17 ), and this is why we started investigating (x-£; + s)I°(x) in the first place! The 

16 evaluators for the "euclidean" product f fgdx, not the "hilbcrtian" f fgdx. 

17 the method was initially developed by the author for the cosine and sine transforms [5, 6] and leads for them 
to the only known "explicit" formulas for C\ for the zero order Hankel transform the problem of computing the 
reproducing kernel had been already solved by de Brangcs [1]. 
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reason is this: the Hubert space of the entire functions .F(s), / <S K a (the Hilbert structure is the 
one from K ai or (jFi, ^2) = ^ Jj{( s ) = i -^"i (s)^2 (s) \r(sj\ 2 ) ver ifi es the de Branges axioms [2], up to 
the change of variable s = \ — iz. Let us recall the axioms of [2] for a (non-zero) Hilbert space of 
entire functions F(z): 

(HI) for each z, evalution at z is a continuous linear form, 



(H2) for each F, z t— > F(z) belongs to the Hilbert space and has the same norm as F, 
(H3) if F(w) = then G(z) = jE^;F( z ) belongs to the space and has the same norm as F 

Let K(z,w) be defined as the evaluator at z: \/F F(z) = (F,K(z, •)). It is anti-analytic in z and 
analytic in w (the scalar product is complex linear in its first entry, and conjugate linear in its 
second entry). It is a reproducing kernel: K(z, w) = (K(z, •), K (w, •)). It is proven in [2] that (HI), 
(H2), (H3) entail the existence of an entire function E(z) with |-E7(z)| > for > 0, such 

that the space is exactly the set of entire functions F(z) such that both and -^|y belong to 
H 2 (3 ! (z) > 0), and the Hilbert space norm of F is ^ J R |F(t)| 2 rg^p - 18 We have incorporated a 2tt 
for easier comparison with our conventions. Then the reproducing kernel is expressed as: 



K M = E{Z)E{W )- E{ C )E{W) d09) 
i(z — w) 

The function E is not unique; if the space has the isometric symmetry F(z) 1— > F(—z), a function 
E exists which is real on the imaginary axis and writing E = A — iB where A and B are real on the 
real axis, the pair (A, B) is unique up to A 1— ► kA, B 1— ► fc _1 _B, A is even and B is odd. If A(Q) 7^ 
(this happens exactly when the space contains at least one element not vanishing at 0) then it may 
be uniquely normalized so that A(0) = 1. Then E is uniquely determined. 

Model examples are the Paley- Wiener spaces of entire functions F(z) of exponential type at 
most t with ||F|| 2 = ^ J R \F(t)\ 2 dt < 00. Then E(z) = e~ iTZ is a possible E function. The 
Paley- Wiener spaces are related to the study of the differential operator — on the positive half- 
line, and an important class of spaces verifying the axioms of [2] is associated with the theory of the 
eigenfunction expansions for Schrodinger operators — j^ + V{x) ([27]). In these examples the spaces 
are indexed by a parameter r (the Schrodinger operator is first studied on a finite interval (0,r)) 
and they are ordered by isometric inclusions (the £?-function of a bigger space may be used in the 
computation of the norm of an element of a smaller space). Typically indeed, de Branges spaces are 
studied included in one fixed space L 2 (R, ^du), are ordered by isometric inclusion and indexed by a 
parameter 19 . Obviously this theory is intimately related with the Weyl-Stonc-Titchmarsh-Kodaira 
theory of the spectral measure. The articles of Dym [14] and Remling [27], the book of Dym and 
McKean [15], will be useful to the interested reader. In the case of the study of 7i we will have 



18 the conditions on F(z) are not formulated in [2] as Hardy space conditions, but they are exactly equivalent. 
19 the axioms allow for "jumps" in the isometric chain of inclusions, as occur in the theory of the Krein strings [15], 
discrete Schrodinger equations being special cases. 
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dv("j) = |r(o + *7)| 2 ^7- It is an important flexibility of the axioms not to be limited to functions 
of finite exponential type, and also the spectral measures are not necessarily such that (1 + 7 2 ) -1 
is integrable. It has turned out in our study of the spaces associated with the 7i-transform that 
the naturally occurring E function is not the one normalized to take value 1 at z = 0. Rather the 
normalization will prove to be lim CT ^ +00 ~^^y^ = 1- This has an important impact on the aspect 
of the differential equations which will govern the deformation of the A' a 's with respect to a: they 
will take the form of a first order linear differential system in canonical form (as generally studied 
in [21, §3].) 

The space of the functions ^(s) = T(s) f(s), f <E K a verify (this is easy) the de Branges axioms, 
with s = | — iz and they were defined 20 in [1]. The spaces K a have the real structure, which is 
manifest in the s variable through the isometry J-(s) i— > J-(s). Rather than with the reproducing 
kernel K(z\,z 2 ) we work mainly with X[s\,S2) — K{—~z\,Z2) which is analytic in both variables. 
Of course then it is X(s, s) which gives the squared norm of the evaluator at s. Writing £ (s) = E{z) 
we obtain from (109): 

X(s 1 ,s 2 ) = (110 

Sl + S 2 — 1 

which is indeed what has appeared on the right hand side of (108). With £(s) = A(s) — iB(s) 1 A 
(resp. B) even (resp. odd) under s i— > 1 — s, this is also: 

^ 1 , S2 )^2-^ ) ^) + - 4 ^)(-^ 2 » (111) 

Sl + S2 — 1 

and for 3?(s) ^ i, < X(s,s) = 2^^^^- so both A and B have all their zeros on the critical 
line. 21 

The method in this chapter has been developed in [5, 8, 6] for the case of the cosine and sine 
transforms, and it leads to the currently only known "explicit" formulae 22 for the structural elements 
£, A, B and reproducing kernels for the spaces for the cosine and sine transforms. So far, almost 
nothing very specific to Ti has been used apart from it being self-adjoint self-reciprocal with an 
entire multiplicative kernel k(xy). The next section is still of a very general validity. 

As was mentioned in the Introduction the realization of the structural elements of the spaces as 
right Mcllin transforms of distributions is a characteristic aspect of the method; the Dirac delta's 
in the expressions for A a (x) and —iB a (x) could have been overlooked if we had only been prepared 
to use functions, and the whole development was based on the computation of (x-rz + s)X s (x) as a 
distribution. This aspect will be further reinforced in the concluding chapter of the paper (section 
9) where it will be seen that the distributions A a (x) and —iB a (x) are very naturally differences of 
boundary values of analytic functions, so they are hyperfunctions [23] in a natural manner. 



20 in the variable z, and associated with the Hankel transform of order zero, rather than with the H transform 

I As x( -s - \£(s)\ 2 ~\e(i- 



21 this is also seen from 2A(s) = S (s) +£(1 - s) as |£(s)| > |f (1 - s)| for K(s) > ±. As X(s, s) - |£(s)1 d£il^£l] 



this is in fact the same argument. 

22 as "explicit" as the Fredholm determinants of the finite Dirichlet kernels are "explicit" 



2G 



Let us consider the behavior of A a (s), B a (s), E a (s) and Tt(E a )(s) for 5R(s) > |. Let us first 
look at E a (s) = ^/a(a~ s + \ f^°(4>t( x ) — 4>a ( x )) x ~ s dx^j . We remark that 4>ti x ) — 0a ( x ) 1S the 
^-transform of -(</>+ (x) + <t>~ (x))l 0<x<a (x). 

Lemma 8. Let k(x) a continuous function on [0, +00) and A G [0,1] be such that ki(x) = 
Jq k(t) dt = 0{x A ) as x — > 00. Let a > and let fix) be an absolutely continuous function on 
[0,a]. Then k(xy)f(y) dy = 0(x A_1 ) as x — > +00. 

There exists C < 00 such that Vx > |fci(x)| < Cx A . Then f£ k(xy)f(y) dy = ±k 1 (xa)f(a) - 
i J a k 1 (xy)f(y) dy, and | /„" k 1 (xy)f'{y) dy\ < Cx A f° y A \f'{y)\ dy. This was easy. . . 

With k(x) = J (2y/x), one has k\{x) = y/xJ\(2-y/x) = 0(xi). We have <j)+(x) — <fia( x ) = 
— J Q a Jo(2y / xy)(0+(a:) + <j>~ (x)) dy and from Lemma 8 this is 0{x~i ). So the integral in the expres- 
sion for E a (s) is absolutely convergent for 3?(s) > |. In particular i£ a is bounded on the critical 
line. But then 7i(E a )(s) = x{s)E a (l — s) is also bounded. Hence: 

Proposition 9. The functions A a and B a are bounded on the critical line. 
Let us turn to the situation regarding 5R(s) = a — > +00. 

Let fix) be a function of class C 2 on [0, a] and e(x) = Jo(2y/xy)f(y) dy. It is 0(x~i). 
There holds ^ae(a:) = ^{^- y yM2^xy))f{y) dy = af{a)J {2^) - f° J (2^xy)yf>(y)dy. Let 
k(x) = Jq Jo(2y / xy) yf'iy) dy. By the Lemma 8 it is 0(x~i). For 5R(s) > |, with absolutely 
convergent integrals: 



af(a) / Jo(2vaa;)a; 3 dx — k(x)x s dx = —ae(a)a s + s eix)x 3 dx (H2) 

J a J a J a 

We show that the left hand side of (112) is 0(a~ s ^) for 5R(s) > |. We apply to k what we did for 
e, £xk(x) = a 2 f{a)J {2^) - / ° M2^y) y(yf')'(y) dy. This is 0(1) (using |J | < 1). So for 
3?(s) > 1, we can compute (-^xk(x))x~ s dx by integration by parts, this gives —ak(a)a^ !S + 
s k(x)x~ s dx. So for Sft(s) > 1 + e we have J°° k{x)x~ s dx = 0(a~ s i). Then regarding 
Ia° Jo(2\fax)x ~ s dx we note that ■j-xJo(2-s/ax) = Jo(2i/ax) — y/axJi{2y/ax), so for 5ft (s) > | + ewc 
can apply the same method of integration by parts, and prove that / a °° Jo{2sfax)x~ s dx = 0(a~ s i). 
So the left hand side of (112) is indeed 0(a~ s ^) for 3?(s) > | and we have: 

Lemma 10. Lei /(#) fee a function of class C 2 on [0, a] and let e(x) = Jo(2^/xy)f(y) dy. One 
has 

e(x)x- s dx = —iae(a) + 0(-)) (8?(s) > J) (113) 
s s 2 

Let us return to Jo(2v / O2?) a;_s = ^(-A)(2a)a 1_s + J^°( Jo(2 v / ax) — v / axJi(2 v / ax))a;^ s dx). 
We want to iterate so we also need x^^/alcJ\{2-Jax) = j2yfax d2 ^ ax 2^/axJi (2^/ax) = ax J${2ifax). 
So we can integrate by parts and obtain that the last Mellin integral is 0(a~ s i) for 3?(s) > | + e. 
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So, certainly: 

{ Jo(2^)x- s dx=—(aJ (2a) + 0(-)) (114) 
J a s s 2 

Using <fi£ = J$ — 7~LPa<f ) t an, l 0a = + ^Pa<i>a arL< l combining (113) and (114) wc obtain: 

Proposition 11. One has for 5R(s) > ctq (Tiere o"o = § /or example): 



a(j) + (a) -Q0 (a) 1 _ a<^(a) 1 

2^ V" A «( s ) = — a ( 1 + — " — + ^^ 

(115a) 

2s + °(^)) -^a(s) = — a (1 __ + (?(_)) 

(115b) 



Theorem 12. One ftas 

lim =f^ = l and ~_ +oc (116) 

So the functions and B a are not normalized as is usually done in [2] which is to impose (when 
possible) to the E function to have value 1 at the origin (which for us is s — 1; the exact value of 
A a {\) will be obtained later.) This difference in normalization is related to the realization of the 
differential equations governing the deformation of the spaces K a as a first order differential system 
in "canonical" form, as in the classical spectral theory of linear differential equations ([21, 10].) This 
allows to realize the self-reciprocal scale reversing operator as a scattering [6]. 

6 Fredholm determinants, the first order differential system, 
and scattering 

Let us return to the defining equations for the entire functions 4>^ and (f>~ : 

4>t +HP a ct>t = Jg (H7a) 
cj>--HP a ct>- = Jg (117b) 

Either we read these equations as identities on (0, oo), or wc decide that TtP a <f>a m f ac t stands for 
Jq Jo(2y/xy)(j)a (y) dy, and the equation holds for i£C; the latter option slightly conflicts with our 
earlier definition of Ti. as an operator on functions or distributions. But whatever choice is made 
this has no impact on what comes next. We shall apply to the equations the operators a-j^ and 
x-§^. As Jg(x) = J (2Vai) we have a-^Jg = xJ^Jg. We write S x = x-§^ + \ = -§^x - \. First we 
have: 

a d~a^ + nPaa d~a^ = J ° + a da~ J ° (118&) 

a da~^ ~ HPaa da~^ = +a ^» (°) J + a Q^ J o ( 118b ) 
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Then, as x£H = -H&x, 5 X H = -HS X , S x P a f = (P a 5 x f) - af(a)S a (x), S x Jg = a£jg + \ Jg: 

S x c/>t ~ HP a 8 x <t>i = (i - a0+(a)) J a + a-^Jg (119a) 

1 3 

Ma + HP a 5 x <t>- = (- + <-(a))J » + a— Jg (119b) 

Combining we obtain: 

(5 9 1 

- + nP a (a—<j>+ ~ M« ) = (<*) + (a) + (120a) 

9 9 1 

a^0a - - HP a (a—<f>- - 5 x cf>+) = +(a#(a) + a0"(a) - -)J Q (120b) 

Comparing with (117a) and (117b), and as there is uniqueness: 

~ = ~ (a0 « (a) + a0 - (a) + ^« (121a) 

The quantity a0+(a) + a<p~(a) will play a fundamental role and we shall denote it by /x(a). 23 So: 

(0^ + ^ + M«M=Ma ( 122a ) 

(o^ + ^ - = (122b) 

It follows easily from this that a^(0+0") = -0+^ + I^CO^) 2 + (#T )*)■ So 

tiMfiW dX = a ^(°)^ W - ^ <^)</>a (*) + ^(0+(a) 2 + C(«) 2 ) (123) 

a i°i ^( a: )C(^)^ = ^(a) 2 (124) 



We then compute: 



cj>+{x)4>-{x)dx= / ((l-D a )- 1 JS)(x)JS(x)dx, (125) 
o Jo 



where we recall 4>+ = (1 + Ha)" 1 Jg , 0~ = (1 - Ha)- 1 Jg, D a = H 2 a . The operator D a acts 
on L 2 (0,a;dx) with kernel D a (x, z) = J Q a Jo(2^/xy) Jo{2^/yz) dy. After the change of variables 
x = at, y = au, z = av this becomes the operator d a on ^(0,1; dt) with kernel d a (t, v) = 
J aJo(2a\^tu) aJo(2ay^uv) du. We compute the derivative with respect to a: 

d f 1 

— / a Jo(2aVtu) aJo(2a\/uv) du (126a) 

f 1 d f 1 d 

= / ((2u— + l)J (2a\^i))aJ a (2ay/uv)du+ / a J Q (2aVtu))((2— u - 1) J (2a\fwu)) du 
Jo ou J du 

(126b) 

= 2aJ (2aVi)J (2aVI>) (126c) 



23 maybe it would be unfair to hide the fact that fJ,(a) = 2a, in this study of H\ In a later section a further mu 
function, associated with a variant of H, will also be found explicitely and it will be quite more complicated. 
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So £d a is a rank one operator, with range CJ (2ay/l)l o<t<i{t). We now use the well-known formula 

logdet(l - da) = -Tr((l - da)- 1 ^) (127) 
da da 

The rank one operator (1 — d a )~ 1 -^d a has the function (1 — d a )~ 1 2a Jo(2ay/i) as eigenvector 

and the eigenvalue is J* Jo(2a^/i)((l — d a )~ 1 2aJo(2ay/v))(t) dt. Going back to (0, a) we obtain 

2 f£ Jo(2 v / aa;)((l — £> a ) _1 Jo(2 v / az))(a;) da; and in view of (125) we have proven: 



d_ 

da 



logdet(l -£>„) = -2 / <!>+ 
Jo 



(x)<t>-(x)dx (128) 



Then, using (124), we have the important formula: 



^{af = -a4-a4- logdet(l - D a ) (129) 
da da 

We shall now also relate <^>+ (a) and <fia( a ) to Fredholm determinants. In fact the following holds: 

a^+(a) = +a4- logdet(l + H a ) (130a) 
da 

acj)-(a) = -a.4- logdct(l - H a ) (130b) 
da 

This is the application of a well-known general theorem, for any continuous kernel k(x, y): if w(x) + 
Jq k(x,y)w(y) dy = k(x,a) for < x < a then w(a) = +^ logdet( . a )(<5(a; — y) + k(x,y)). A proof 
may be given which is of a somewhat similar kind as the one given above for (128), or one may 
more directly use the Fredholm's formulas for the determinant and the resolvent. 2 The theorem 
is proven in the book of P. Lax [20], Theorem 12 of Chapter 24 (Lax treats the case of a kernel on 
(a, +oo), here we have the simpler case of a finite interval (0, a).) This means that /x(a) has another 
expression in terms of Fredholm determinants: 

d det(l + Hg) 
^ a) = a d- a l ° g det(l-H a ) (131) 

Combining (129) and (131) we obtain: 

„ d d, , , . / d, dct(l + ff a )\ 2 d d, det(l + H a ) ,., . 

- 2a Ta a Ta + g °> = [ a d-a ^ det(l - h]) ) ~ Vda ^ ^lljj ^ 

2a d a d logdet(l H ) - ( a - log det ^ + ^ V + a d a - log det ^ + ^ (132b) 
da da a \ da det(l — H a ) J da da det(l — H a ) 

2a4 L acf>+(a) = -fJ,(a) 2 + a// (a) (132c) 
da 

2a4 L a6~ (a) = +/i(a) 2 + ou'(a) (132d) 
da 

^a(0-(a) - 0+(a)) = a(0+(a) + 0"(a)) a (132c) 

These Fredholm determinants identities arc reminiscent of certain well-known Gaudin identities [22, 
App. A16], which apply to the even and odd parts of an additive (Toeplitz) convolution kernel on 



24 let us recall that for a continuous kernel on a finite interval, the formula of Fredholm for a determinant as a 
convergent series always applies, even if the operator given by the kernel is not trace class, which may happen. 
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an interval (—a, a); here the situation is with kernels k(xy) which have a multiplicative look, and 
reduction to the additive case would give g(t + u) type kernels on semi-infinite intervals. 

We have defined A a = ^0+ + H(f>t) and -iB a = ^{-(f>~ + H<p~). Let us recall that 
here </>J is restricted to [0, +00) and is then tempered as a distribution. Using the differential 
equations (122a) and (122b) and the commutation property S x Ti = —7i.5 x , 5 X = x-j^ + i, we 
have 5 x A a = ^(5 X 0+ - H8 X 4>+) = ^(a£ + \ - p(a))(0- - W#T) = -(a£ - fi(a))(-iB a ) and 
S x (-iB a ) = ^(-(o& + § + ii(a))<l>+ - (a& + § + K^Wt) = -(a& + »{a))A a . The following 
first order system of differential equations therefore holds: 

d 

a—A a = -fi,(a)A a - 5 x {~iB a ) (133a) 
aa 

a^-{-iB a ) = +fi(a)(-iB a ) - 5 x A a (133b) 
aa 

Then we also have the second order differential equations (aj^ ~ L i { a )){ a ~§^ +A*( a ))^a = +^ x A a and 
( a 3a f J, ( a ))( a "5a — M( a ))( — *-^a) = +5 x (—iB a ) 1 or, taking the right Mellin transforms, and writing 
s = \ + i*y, S x = ij: 

-a—a—{-iB a ) + {^{af + a/V (a)){-iB a ) = 7 2 H^) (134b) 

With the new variable u = log(a) we obtain Dirac and Schrodingcr equations which are associated 
with this study of TL, modeled on the study of the cosine and sine transforms summarized in [5, 6]. 
All quantities in the statement of the theorem will be completely explicited later in terms of Bessel 
functions, but we keep the notation sufficiently general to allow, if an interesting other case arises, 
to write down the identical results: 

Theorem 13. For each a > let and <p~ be the entire functions which are the solutions to: 

d>t(x) + [ Jo{2^xy)^(y) dy = J (2^) (135a) 
Jo 

<f>-(x) - f J o (2V^y)0a (V) dy = J (2v^) (135b) 
Jo 

Let H a be the integral operator on L 2 (0,a]dx) with kernel Jo(2^/xy). There holds: 

^+(a) = +-^-logdet(l + ff a ) (135c) 
da 

^-(o) = -^-logdet(l-fr„). (135d) 
da 

The tempered distributions A a = ^(1+W)(0+ lo<x<oo) and B a = i^-(—l+Tt)((f)~lQ <x<00 ) vanish 
on (—00, a) and are respectively self -reciprocal and skew-reciprocal under H.. Their completed right 
Mellin transforms A a (s) = T(s)A a (s) andB a {s) = T(s)B a (s) are entire functions with all their zeros 
on the critical line, they are respectively even and odd for s <-> 1 — s, and they verify the following 
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Dirac and Schrddinger types of differential equations in the variable u = log(a), — oo < u < +00, 



d 

au 


= -n(a)A 


±B 

au 


= +v(a)B c 




= (-- 

\ du 2 


2 n 

1 B a 


( d 2 
~ V du 2 


V + {\oga) 


= M(a) 2 " 


V-{\oga) 


= K a ) 2 + 


fj,(a) 


d lo d 

du ^ d 



duja) = _ 2 d 2 logdet(l + H a ) 

du 
dfi(a) 



du 2 

d 2 l0gdct(l -H a ) 



du 



du 2 

- a</>+ (a) + a</>~ (a) 



(135e) 
(135f) 
(135g) 

(135h) 
(135i) 
(135j) 
(135k) 



where s 



17. 



Let us consider the Hilbert space of pairs 



a(u) 
(3{u) 



on R with squared norms \a(u)\ 2 + 



(u)| 2 4r, and the two equivalent differential systems in canonical forms: 





' 1" 


d 


u(e u ) 


( 


-1 


du 


u{e u ) 





a(u) 


= 7 


a(u) 


) 


p(u)_ 


p{u)_ 



( 


' 1" 


d 


~-fi(e u ) 




-1 


du 


fi(e u ) 



a(u)+/3{u) 
-a(u)+j3(u) 



7 



a(u)+/3(u) ' 
-a(u) + (3(u) 



The components obey the corresponding Schrddinger equations: 

-a"(u) + V+(u)a(u) = 7 2 a V+(u) = u{e u f - 

-p"(u) + V-{u)(3{u) = 7 2 /3 V-{u) = u(e u ) 2 + 



du{e u ) 

du 
dfi(e u ) 

du 



(136) 
(137) 

(138a) 
(138b) 



Regarding the behavior at —00, we are in the limit-point case for each of the Schrddinger equations 
(138a) and (138b) because clearly (say, from the defining integral equations for <^>+ and <fi~) one has 
'fit ( a ) ~~ *o— >o ^o(O) = 1, 4>a (a) —f a ^o 1, fJ-(a) ~ a — >o 2a, so the potentials are exponentially vanishing 
as u — > — oo. Perhaps we should reveal that one has exactly u{a) = 2a = 2e u so we are dealing here 
with quite concrete Schrddinger equations and Dirac systems whose exact solutions will later be 
written cxplicitely in terms of modified Bessel functions, but we delay using any information which 
would be too specific of the 7i-transform. 



For each 7 e 



«4exp(u)(| +*7) 
.^exp(«)(5 +«7) 



(139) 



is a (non-zero) solution of the system (136), and we now show that it is square-integrablc (with re- 
spect to du — dlog(a)) at +00. Let us recall the equation (111) (s+z — l)X a (s, z) = — 2iB a (s)A a (z)~ 
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2iA a {s)B a (z), from which we deduce 

a-^-X a (s,z) -- 
oa 



-2A a (s)A a (z) - 2(iB a (s))(iB a (z)) 



(140) 



We have 25 \\X?\\ 2 = X a (s,s), A a (s) = A a (s), iB a {s) = iB a (s), so 



a-^\\X*\\ 2 = -2\A a (s)\ 2 -2\B a (s) 



(141) 



and as of course lim a ^ +oc HA""!) 2 = (we have ||<-f s a || 2 < \X^\ 2 (x) dx for a > 1) we obtain: 



Vs e 



\X° 



\Ba( S )\ 2 )- 

a 



(142) 



This establishes the square-integrability at +oo of 



, for any s £ 
and 



The solutions of (136) with eigenvalue 7 = are • Aa ^2> and ^(i)- 1 ■ The former is square- 
intcgrablc, so from 2 < t + t^ 1 the latter then necessarily is not. This confirms that the Dirac 
system (136) is in the limit point case at +oo (according to a general theorem of Levitan [21, §13, 



a(u) b(u) 
c(u) d(u) 



with continuous coefficients is 



Thm 7.11 any first order differential operator [ J 1 -jr- 
in the limit point case at infinity). So the pair (139) is in fact, for any 7 £ C, the unique solution of 
(136) which is square-integrable at +00. Also the Schrodinger equation (138b) is in the limit point 
case as not all of its solutions are square integrable at +00. Whether the limit-point case at +00 
holds for equation (138a) is less evident. Let us recall from [10, §9, Thm 2.4] and [26, §X, Thm X.8] 
that a sufficient condition for this is the existence of a lowerbound liminf M _> +00 V+(u)/u 2 > —00. 
We will prove in the next chapter that /i(a) = 2a = 2e u so this is certainly the case here. In the 
present chapter only the fact that the Dirac system is known to be in the limit-point case will be 
used. 

We now take uo = logao and apply on (u , 00) the Weyl-Stone-Titchmarsh-Kodaira theory ([10, 
§9], [21, §3]). Let tp(u, s) be the unique solution of the system (136) for the eigenvalue 7, s 



and with the initial condition ipi u o, s) 



condition 0(ito,s) = 



- roi 



+ 17, 

and let 4>(u, s) be the unique solution with the initial 
Let m(7) - 0(M, s) + c/)(u, s) for $57 > be the unique solution which is 
square-integrable on (uq, +00). So 771(7) = b"°(s) ' s = 3 +*7> ^( s ) < \- It is a fundamental general 
property of the m function from Hermann Weyl's theory that 3(771(7)) > (for 3(7) > 0.) Here, 
we have a case where the m-function is found to be meromorphic on all of C; so we see that its 
poles and zeros on R are simple. Furthermore, the spectral measure v is obtained via the formula 
via, b) = lim e _ > o+ — / 3m(7 + ie) 07 (under the condition v{a, 6} = 0). We obtain: 

-<5(7 - 3p) 



(143) 



The spectrum is thus purely discrete and the general theory tells us further that the finite linear 



combinations ^ 



tp(u,p) have squared norms 



P -iB' ao (p)^P 



c p \ 2 and also that they are 



25 let us recall the notation Xf = T(s)X% S L 2 (a, +00; dx). 



33 



dense in L 2 ((u , oo) -> C 2 ; du). For B ao (p) = 0, ip(u,p) = A ao (p) 1 k^"'^] 1«>« (u), so the 



vectors 



2Axp (u) (p) 
2B eX p (u )(p) 



cxp(u) (p) _ 

u>m ( u ) are an orthogonal basis of L 2 ((uo,oo) — > C 2 ; ^c?u) and they 
satisfy ||^p°|| 2 = — 2.4 ao ( ( o) iB' ao (p). Similarly a spectral interpretation is given to the zeros of A ao 
if one looks at the initial condition [ J ] . The factors of 2 and \ , have been incorporated so that the 
statement may be translated (taking into account results established later) into the fact that the 
evaluators K ao (p, z), for B ao (p) = 0, are an orthogonal basis of the Hilbert space of the functions 
r(z)/(z), / G K aa . This last statement is a general theorem (under a certain condition) for spaces 
with the dc Brangcs axioms [2. §22]. 

To discuss in a self-contained manner the generalized Parseval identity which is associated with 
the differential system on the full line, it is convenient to make a preliminary majoration of ||A""|| 2 , 



sft(s) = i. From (108) we have, for Ofc(s) = \: \\X%f = 2M(£ a (s)E' a (s)). And £ a (s) = T(s)E a (s). 
And E a (s) = Y / a^a _s + 5 f^i&tix) - ^ ( x ))x~ s dx\ . We know from the discussion of Lemma 9 that 
the integral in the expression for E a (s) is absolutely convergent for 5R(s) > j. Hence by the Ricmann- 
Lebesgue lemma E a (^ + i'y) ~ a~ n as j — > oo, 7 € R. Similarly, E a (| + ij) ~ — log(a)a _i7 . So, 
with ll^fH 2 = |r(s)| 2 || A""|j 2 and using Stirling's formula we obtain: 

Lemma 14. For each given a > one has ||A"|| 2 ~ 2 log |s| as \s\ — ► oo, 5R(s) = |. 

From (104) expressed using A a and B a we see that Ba ^f > is square integrable, so s~ 1 B a (s) is 

s 2 

square integrable on the critical line (with respect to \ds\). Then using again (104) we see that 
s~ 1 A a (s) is also square integrable on the critical line. 26 Let us pick a function F(s) on the critical 
line which is such that sF(s) is square integrable. Then F(s)A a (s) and F(s)B a (s) are absolutely 
integrable on the critical line and (f^^—i \F{ s )A a {s)\ ^r) 2 < C Jjf( s ) = i ^ an( l similarly 

with B a . If we define a F (u) = 2 / R(a)=i F(s)A a (s)^ and /3 F (u) = 2 / R(a ) = £ F(s)B a (s)^- we 
then compute: 

r(K( S )i 2 + is;( S )i 2 )^ | d 



|ai,(w)| 2 + |/3 F (M)| 2 du<C 

«o JSR(s) 



£1 

2tt 

all 2 



"/ !M <co (144) 



So ai?(w) and Pf(u) are square integrable at +00. More precisely the above upper bound holds as 
/«(.)=* I^WAiWI^ and 



well for L w=i |.F(s)4,(*)|lf2L and L, . =i |F(a)S (s)|i^L. So the double integrals 



//„ 



Axp(«)(2)A X p(n)(s)-? : '(s) JwL 2 du ( 145a ) 

ti <"<oo,!R(s) = ± /7r lM s JI 

S exp( „)(2;)B ex p(„)(s)J"(s) |f/L 2 ^ (145b) 

M <u<oo,SR(s)=i ' !7r l i l s ;l 



'we know in fact according to proposition 11 that A a and _B a arc bounded on the critical line. 
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where z G C is arbitrary, and T(s) = T(s)F(s), are absolutely convergent and Fubini may be 
employed. Using (140): 



X exp ( Uo )(z,s) = 2 / (A cxp ( u ){z)A cxp ( u ){s) + B exp ( u )(z)B CX p(u)(s))du (146) 



And we obtain the following identity of absolutely convergent integrals, for any J-(s) = T(s)F(s) 
with S F(,s)GL 2 (5R(.s) = i;^i): 



^oxp(u )( 2; : S )- ?r ( S ) 27r |r(l)| 2 = 



We shall prove that this identity holds under the weaker hypothesis F(s) G L 2 ($i(s) = \; -^-). 
First, still with sF(s) square integrable we suppose additionally that F = f with / G -Kexp(uo) 27 - 
The hilbertian kernel K cxp{uo) (z, s) is X cxp{uo) (z, s) so K CKp{uo) (z, s) = Af exp(uo) (z, s). The equations 
give then: 



T{z) = 


y ( 




Ju a 


ap(u) = 


2 I 




Jnt(. 


Pf{u) = 


2 I 







(148a) 



WW o.L'm2 ( 148b ) 



^)ga(«) ^T7w^ia ( 148c ) 



H 

'27r|r(s) 

Jds| 
'27r|r(s) 



We have worked under the hypothesis that sF(s) is square integrable. To show that the formulae 
extend in the L 2 sense, we first examine: 

Mu)\ 2 = 4 / W( Sl ) 9 |ff [ / ^)A( g2 ) 9 |ff' (149) 

y«(«0=i 27r l r (si)l Jg?( S2 ) = i 2tt\T{s 2 )\ z 

\a F (u)\ 2 + \(3 F (u)\ 2 ^= II ^i)^)A^, (uo) ( ai ,g a ) Iff 1 If 881 (150) 

2 JJu( Sl )=^ 2tt |r(si)| 2 2tt |r(s 2 )| 2 

There was absolute convergence in the triple integral used as an intermediate. Also X exp i UQ \(si 1^2) = 



Hence: 



X oxp{ua ){s 2 ,s 1 ) and J u(si)= i F(si)X exp ( Uo )( s 2> g i) 27 r |r(s' )|^ = -^faa) 

r(M«)i 2 +^(«)i a )^= / if( S )i 2 ^= r i/wi 2 ^ (i5i) 

So with an arbitrary / G K ai F = f, J-'(s) = r(s)/(s), the assignment / 1— > (cufiPf) exists in 
the sense of L 2 convergence when one approximates / by a sequence f n in K a such that sf n (s) is 
in i 2 (5R(s) = |; ^p-), and / 1— > (o.f,Pf) is linear and isometric. We check that its range is all of 
L 2 (ito,oo; 4r) © i 2 (uo:Oo; 4r). For this let us identify the functions a w (u) and (3 w (u) which will 
correspond to !F(s) = X ao (w, s) (ao = exp(uo).) On one hand from (147) it must be the case that: 



\ds[ 



Vz G C / ^ ^(2:, 5)^(^,5) ^ ]r/ ^ }2 = / (Ax P («)(2 : )a«.(u)+Soxp(t l )(z)/3, 1 ,(w))du (152) 



27 this is certainly possible as we know that the f(x) which are smooth, vanishing on (0, a) and of Schwartz decrease 
as x — > +00 are dense in K a . 
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The left hand side is X ag (w, z) which on the other hand is given by the formula 2 (A a (z)A a ('w) — 
B a (z)B a (w)) du. The functions u i— ► ' z ^ ^ are cer tainly dense in L 2 ((uo,oo) — > C 2 ; 4^) as 

we know in particular that the pairs for the p's such that B ao (p) = give an orthogonal basis. So 
we have the identification on (ito, +00): 

a w (u) = 2 Axp(u) (w) Pw{u) = -2B exp („)(w) (153) 

This proves that the range is all of L 2 ((uo,oo) — > C 2 ; Let us note that in this identifica- 

tion the hilbcrtian evaluator K ao (w, •) is sent to the pair u > 2 l«> Uo ( u )(^a(w) , B a (w)). To 
check if all is coherent we compute the hilbertian scalar product (K ao (w, •), if ao (z, •)). We obtain 
4 -C AHAO) + £ a (»£ Q (z)4f = 2 /~ - B a (w)B a (z) du = X exp(uo) (w, z), which is 

indeed K ao (w, z). 

Let us return to the consideration of a general F(s) 6 L 2 (5R(s) = i; ^-). Under the hypothesis 
that s F(s) is square integrable we have assigned to F the functions 

M«) = 2 / ^)A( S ) 9 , 2 = / F(s)2A a (s) 1 -^ (154a) 

/»(«)=! 2tt |r(s)| 2 Jsr (s)= i 2tt 



M«) = 2 / ^( s ) Sa ( s )_L_L_^ = / F(5) 2B ( S ) L_L (154b) 



which arc square-integrable at H-oo. From (147) there holds, for any ao = cxp(izo): 

X exp{uo) (z,s)F(s) l —^= / (2A exp(u) (z)a i r(u)+2B^; ) (z)^ F (M))— (155) 
5R(s)=i 27T J„ 2 

The function of z on the left side is the orthogonal projection F ao of F to the space K ao . So, 
we deduce by unicity of (w)l«> tl0 (u) = Q ? i? UQ (w) and /3f(u)l tl > Uo (it) = /3f u (u). We then obtain 
ll-^a || 2 = XT°(l a F(M)| 2 + \Pf(u)\ 2 )^- so q;f and /3f are square-integrable on (—00, +00), and as 
U a K a is dense in L 2 (0,oo;dx) the assignment F t— > (oif,/3f) is isometric, and also it has a dense 
range in L 2 (M — > C 2 ; We can then remove the hypothesis that sF(s) is square integrable and 
define the functions oif and (3f to be the limit in the L 2 sense of functions a n and (3 n associated 
with FnS such that \\F — F n \\ — *■ and the s F n are square-integrable. Summing up: 



Theorem 15. There are unitary identifications L 2 (0, 00; dx)— ^X 2 (5R(s) = ^; ^-)^L' 2 



-'2. du \ 

' > 2 > 



given in the L? sense by the formulas, where ^(s^ - 



2 ■ 



= /(*) = / f(x)x~ S dx f(x) = / F^x*' 1 ^ (156a) 

JO J5R(s) = I 

a(«)= lim / F n ( S )2i^; ) ( S ) ^ (F„ ^ F; S F„(.s)eL 2 ) (156b) 

£(«)= hm / F n (s)2C^)W^ (156c) 

F(*)=lim/ (a(u)2 J 4 exp(u) (s)+/3(u)2S exp(u) (s))— (156d) 



2- 
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The orthogonal projection of f to K ao corresponds to the replacement of a(u) by a(u)l u>Uo (u) and 
of (3{u) by P(u)l u>UQ (u) (u = \og(a a ).). The unitary operators f i-> H(f), F(s) i-> x(s)F(l - s), 
correspond to (a, (3) i— > (a,—f3). For f = one has a(u) = 2j4 e3C p( u )(e)l u> i O g( O0 )(ii) and (3(u) = 
— 2B oxp ( tl ) (z)lti>iog(a ) ( u )- The self-adjoint operator F(s) i— > 7-F(s) (s = | + *7^ corresponds to the 
canonical operator: 

^(e* 


















^(e u ) 



(156e) 



which, in L 2 (R — ► C 2 ; 4^), is essentially self-adjoint when defined on the domain of the functions of 
class C 1 (or even C°° ) with compact support. The unitary operator e lrH acts on L 2 (0,oo;dx) as 
f(x) i-> e^ T f(e T x). 

For the statement of self-adjointness we start with a and j3 of class C 1 with compact support, 
define F by (156d) and integrate by parts to confirm that r yF(s) corresponds to H {[%]). We know 
by Hermann Weyl's theorem that in the limit point case the pairs [p] of class C 1 with compact 
support are a core of self-adjointness (c/. [21, §13].) On the other hand we know that multiplication 
by 7 on i 2 (5R(s) = ^; -^) with maximal domain is a self-adjoint operator. So the two self-adjoint 
operators are the same. 

Having discussed the matter from the point of view of the isometric expansion we now turn to 
another topic, the topic of the scattering, or rather total reflection against the potential barrier 
at +oo. Another pair of solutions of the first order system (136) (hence also of the second order 
differential equations) is known. Let us recall from equations (95a), (95b) that we defined j a = 
\/a{5 a - (j)+\ <x<a) = \faH<t>t and ~ ik a = \fa{8 a + 4>a 1 o< K <a) = y/aH(j)~ . Using again (122a) and 
(122b) it is checked that j a and k a verify the exact same differential system as A a and B a : 

d 

a^rja = -n(a)j a + i5 x k a (157a) 
oa 

d 

a—k a = +n(a)k a - i8 x j a (157b) 
oa 

The right Mellin transforms j a (s) and —ik a (s) are defined as 

^ pa 

j a {s) = a?~ s - y/a (j)+(x)x~ s dx (158a) 
Jo 

ra 

-ik a (s) = a^ s + y/a (j)~{x)x~ s dx (158b) 
Jo 

As 4)^ and 4>~ are analytic these are meromorphic functions in C with possible 28 pole locations at 
s = 1, 2, ... From the point of view of the Schrodinger equation (138a) and as u = log(a) — > — oo, 
a — ► 0, we thus see that, for s = ^ + ?7, 7 € K, ia(i+*7) and j a (h — vy) are two (linearly independent 
for 7 ^ 0) solutions, differing from e~ tJU and e' 7 " by an exponentially small (in u = log(a)) quantity 
(and similarly with — ik a with respect to the Schrodinger equation (135h)). So we have identified 
the unique solutions which verify the Jost conditions at — oo. 



'the poles do exist. 
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As P a 4>a is square integrable, and also using (85a), we have on the critical line: 

r(*)£(s) + T(l - s)f a {\ - s) (159a) 



= r (s).7a(s) - 


1-1(1 — sja 


= mla(s) - 


|-r(l - s)a s 




(- T(l - s)a s 


= r(.s)a^- s - 


r- T(l - s)a s 



(159b) 
(159c) 
(159d) 



+ r(s)V^/ (0+(.t) - J (2V^))x- s da; (159e) 



As Jo(2y / ox) — 4>a(x) is square integrable both integrals are simultaneously absolutely convergent 
at least for \ < 5R(s) < 1 (the | can be improved, but this does not matter). As the boundary 
values on the critical line coincide we have an identity of analytic functions. We recognize in 
J°° J$(2ijax)x~ s dx, which is absolutely convergent for 3?(s) > |, the quantity g s (a) (equation 
(73)). And from equation (74) we know g s {a) = x(s)a s_1 — Jo(2yfax)x~ s dx. So 

/>oo 

T(s)T a (s) + r(l - - s) = T(s)ai- S + r(a)V5 / <tf(x)x" s , (160) 



which is indeed 2.A a (s). From the equation (158a) j a (s) = a^(a~ s — <f)+(x)x~ s dx) (valid as 
is for 5R(s) < 1) the function u i— * j Q (s) differs from it i— > e~ t71i by an error which is relatively 
exponentially smaller (we write s = -| + vy, ^(7) > —5). So j is the Jost solution at — 00 of the 
Schrodinger equation (135g). The identity relating B a (s) and k a {s) = ia?(a~ s + f£ 4>~{x)x~ s dx) 
is proven similarly. 

Theorem 16. The unique 29 solution, square integrable at u = +00, of the Schrodinger equation 
(135g) (resp. (135h); 7 7^ 0) is expressed in terms of the functions j a {\ *t) ( res P- ~ik a (^+ij)) 
satisfying at —00 the Jost condition j a ~ u _^_ oa e ~ % i u (resp. —ik a ~ u _^_ oa e~ nw ) as: 

A a (s) = \ (T(s)f a (s) + T(l - s)3~ a (l - s)) (161a) 

s«(«) = \ (r(s)fcT(s) - r(i - s )T a (i - «)) (ieib) 

Let us add a time parameter t and consider the wave equation: 

92 V-£W«)=0 (162) 



dt 2 du 2 ^ du r 

Then $(i,u) = e l7 *j exp („) (| +17) is a solution which behaves as e %1<yt ~ u ^ asw^ —00. This wave is 
thus right-moving, it is an incoming wave from u = — 00 at t — —00. For a given frequency 7 there 



29 here we make use of the fact that (135g) is in the limit point case at +00, because it is proven in the next chapter, 
or known from (66a), (66b), that fi(a) = 2a = 2e u , in this study of the Ti. transform. 
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is a unique, up to multiplicative factor, wave which respects the condition of being at each time 
square integrablc at u = +00. This wave is e l7 *„4 oxp ( u ) (| + 27). So the equation (161a) represents 
the decomposition in incoming and reflected components. There is in the reflected component a 
phase shift B 1 = arg%(s), the solution behaving approximativcly at u — > —00 as C(j) cos(ju + \0^). 
This is an absolute scattering, as there is nothing a priori to compare it too. We will thus declare 
that equation (162) has realized x(s) as an (absolute) scattering. Similarly the Schrodinger equation 
(135h) realizes — x( s ) as an absolute scattering. 

We have 2A a {\) = 2Y{\%{\) and j a {\) = l-o* / a #(a:)ar* dx. So lim^ M\) = = 
y/ir. On the other hand a-£A a (^) = -n(a)A a {\) and fj,(a) = a-£ log ^(1-^°) - so: 

, /In /—^r~As 1— dct(l — H a ) 

^2) = ^2) = ^d^rr4 (163) 

We have a^||A?|| 2 = -2A a (\f '■ And Af| = Y{\)X\. So: 

Theorem 17. The squared-norm of the evaluator f 1— > X'i (/) = on the Hilbert space 

K a of square integrable functions vanishing on (0, a) and with Ti transforms again vanishing on 
(0, a) is: 

where H a is the restriction ofTi to L 2 (0,a]dx). 

It will be seen that det(l + H a ) = e a ~^ a and det(l — H a ) = e~ a ~h a . Having spent a long time 
in the general set-up we now turn to determine cxplicitcly what the functions (f>^ , 4>~ , etc. . . are. 



7 The i^-Bessel function in the theory of the TL transform 

Let us recall that we may define the H transform on all of L 2 (R; dx) through the formula Ti{f ) (A) = 
t/(^)- This anticommutes with f(x) 1— > /(— x), and Ti. leaves separately invariant L 2 (0, +00; dx) 
and L 2 (— 00, 0; dx). We defined the groups r a : f(x) 1— > /(x — a) and r^f = H.T a TC. We observed 
that the two groups are mutually commuting, and that if the leftmost point of the support of / is 
at a(f) > then the leftmost point of the support of rf(f), for any b > 0, more precisely for any 
b > —a(H.(f)), is still exactly at a{f). From this we obtain the exact description of K a : 

Lemma 18. One has K a = T a TfL 2 (0, +00; dx). 

Let now Q be the orthogonal projection L 2 (R; dx) 1— ► L 2 (0, +00; dx). The orthogonal projection 
Q a from L 2 (0,oo;dx) to i4T a is thus exactly T a T#QT- a T^t a . It will be easier to work with R a = 
Qr^ a T^ a , especially as we are interested in scalar products so we can skip the T a r# isometry. First, 
we obtain gf(x) — R a (ft)(x), for ft(x) = e~ tx . The part of r_ a (/ t ) supported in x < will be 
sent by r^ a to a function supported again in x < 0. We can forget about it and we have thus first 
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e~ ta e~ tx l x> o(x), whose Tt transform is e~ ta 7 exp(— |), which we translate to the left, again we cut 
the part in x < 0, and we reapply ft, this gives gf (x) = e~ a ( t+ ^e~* a: l 2 ;>o(a;). In other words we 
have used in this computation: 

Qr-aT* a = HQT- a HQT- a {a > 0) (165) 

The orthogonal projection f" := Q a (ft) of ft(x) = e~~ tx l x>0 {x) to K a is thus T a T#{gf). We can 
then compute exactly the Fourier transform of / t a as ft {it) = {e~ TX , T a rf (<7"))l 2 (r) which is also 
(r # a r_ a e— ,g?) L 2 {U ) = {9r,gf) = e-^+De-^+i)^. Hence: 

Lemma 19. TTie orthogonal projection / t a to if a 0/ e~ tx l x> o{x) has its Fourier transform / ( a (A) 
which is given as: 

ftfir) = e - ait+ T+ T +i) -L- (166) 

t + T 

The Gamma completed right Mellin transform Tt{s) of /" is the left Mellin transform of ft {it). 

poo s — 1 

f?(x)X?(x) alx = F?{s) = e - a(t+ 7) / e - a(T+ ^ dr (167) 

Jo * + T 

Let us write for the element of L 2 {0, +00; alx) such that T a T,f = X* . We have F?{s) = 
{X?, / t a ) = (W s a , .9?) = e - Q ( t+ T) J °° W s °(:r)e- te dz. So the Laplace transform of W£{x) is exactly: 

W s a (x)e- te d:i; = y e "° (T+v) ^: dr ( 168 ) 
Writing 7^ = J Q e -( t + T ) x dx, we recover W"(x) as: 

/•CO 

W s a (x)=/ e-°( r+ - V^e - ™ dr (169) 
Jo 

Then we obtain W"(x)W"(x) dx which is nothing else than X a {s, z): 

Theorem 20. The (analytic) reproducing kernel associated with the space of the completed right 
Mellin transforms of the elements of K a is 

X a {s,z)= ff e -°(«+l+"+£) g^!Z^ dtdu (170) 

JJ[0,+oo) 2 t + U 

Here is a shortened argument: the analytic reproducing kernel X a {s,z) is the completed right 
Mellin transform of X?(x), so this is / °°(# s o , e'^t 8 ' 1 dt. But for U{s) > |, {X?,e~ tx ) = 
T{s){Q a {x~ s l x>a ),e~ tx ) = T{s){x~ 3 l x>a ,f^) = J^(s) {Q a is the orthogonal projection to K a ). 
This gives again (170). 

To proceed further, we compute {s + z — l)X a {s, z). Using integration by parts, multiplication 
by s (resp. z) is converted into — t-i (rcsp. —u-j^; there are no boundary terms.) 

sXJs,z)= [[ {a{t --) + -!—)e- a{t+ i +u+ i ) - — 1 dtdu (171a) 

JJ[0,+oo)' t t + U t + U 

zX a {s,z) = ff ^{u--) + ^)e- a ^ +u+ ^ i —^—dtdu (171b) 

JJ[0,+oo)2 U t + U t + U 
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(s + z- l)X a (s, z) = a[[ (t - \ + u - I) e -(*+i+«+i)^l- 

Mo,+°°) 2 t U t + U 



■ dtdu 



oo />oo 



-a(t+\) t s-l dt / e -a(u+i) u z-l du 
JO 



a(t+l) t s-2 dt / e -a(n+i) uZ -2 du ( mc ) 



The X-Bessel function is K s (x) = \ J °° e~ x ^ t - t+ ^h s ^ 1 dt = J °° e _:rcoshtl cosh(sit) du. It is an even 
function of s. It has, for each x > 0, all its zeros on the imaginary axis, and was used by Polya in a 
famous work on functions inspired by the Riemann ^-function and for which he proved the validity 
of the Riemann hypothesis [24, 25]. We have obtained the formula 

s + z — 1 

To confirm £ a (s) = 2^/EK s {2a), let us define temporarily A(s) = \-Ja J °° e - a(t+ ^ (1 + i)** -1 
and —iB(s) = \^fa / °° e~ a ^* + ^(l — \)t s ~~ x dt which are respectively even and odd under s i— > 1 — s 
and are such that £?(z) = A(z) - iB(z). We have Vs, z <G C -iB(s)A(z) + ^4(s)(-iS(z)) = 
— iB a (s)Aa(z) + A a (s)(—iB a (z)) and considering separately the even and odd parts in z, we find 
that there exists a constant k(a) such that A(s) = k(a)A a (s) and B(s) = k(a)~ 1 B a (s). Let us check 
that linio-^oo ■jjay = 1- It is a corollary to linv^oo K (T (x)/K (T ^.i(x) = which is elementary: 
cxp(— x coshu)e au du — O(l) (a — > +oo), and for each T > 0, exp(— a; coshM)e crtl e?u > 
Tcxp(— xcosh3T)e 2crT , J" T exp(— xcoshu)e au du < Te aT , and combining we get K a {x) = (1 + 
o(l))| exp(— xcoshu)e CT " du. So limsup (T _ >+00 j^"^].) < e~ T for each T > 0. Using (116), we 
then conclude k(a) = 1. 

Let us examine the equality £ (s) = 2yfaK s {2a) = y/a J °° exp(— a(t + i))t s_1 c?t. It exhibits 
£ a as the left Mellin transform of -yaexp(— o(t + j)), so the distribution E a is determined as the 
distribution whose Fourier transform is y/aex~p(i a(X — A -1 )). Using r a and rf , this means exactly: 

E a = y/i t* T a 5 = y^ Hr a HS(x - a) (173) 

We exploit the symmetry K s = K- s , which corresponds to E a (X) = E a (— A -1 ) = —iXHE a {\), so 
the unexpected identity appears: 

E a = ^HE a (174) 



From (173) wereadK-Ea = ■\/aT a Jo(2y/ax) = v / aJo(2-\/a(a; — a))l x>a (x). Using (174), and recalling 
equations (96c) and (96d) we deduce: 

0+(z)+0-(x) = Jo{2 ^ a{x _ a)) = /o(2v / a(a _ ~) (175a) 



J (2^a{x-a)) = —I (2^/a(a~x)) (175b) 



2 dx dx 

We knew already from equations (66a), (66b)! Summing up we have proven: 
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Theorem 21. LefH be the self-reciprocal operator with kernel Jo(2^/xy) on £ 2 (0, oo; dx). Let H a be 
the restriction ofH to L 2 (0,a]dx). The solutions to the integral equations + H a <f>+ = Jq(2^/ox) 
and — H a (/)~ = Jo(2y/ax) are the entire functions: 



d 



<t>t{x) = I {2y/a(a - x)) + —I {2^a{a-x)) 



(176a) 
(176b) 



4>~(x) = I (2y/a(a - x)) - —I {2^a{a-x)) 

One has 1 — a = (a) = ^ logdet(l + H a ) and 1 + a = <j)~ (a) = — ^ logdet(l — H a ), and 

det(l + # a ) = e +a -i a2 (176c) 
det(l - H a ) = e- a -^ (176d) 

The tempered distributions A a = -^(1 + Ti)<f>a an d —iB a = 2 ( — 1 +^)0a > respectively invariant 
and anti-invariant under Ti, are also given as: 

d 



Mx) = *y (S a (x) + l x>a (x)(j (2y/a(x-a)) + — J {2^a{x-a)))^ 



-iB a {x) = ^ (S a {x) - l x>a {x)(j (2^a(x-a)) - ^J {2^/a(x - a)))) 



(176e) 
(176f) 



Their Fourier transforms are L e lA:E ^4 a (a;) dx = -^(1 + jr) exp(ia(A— j-)) and —i J R e l * x B a (x) dx 
■^(1 — 4-)exp(ia(A — 4)). 27ie Gamma completed right Mellin transforms are: 



T(s)T a (s) = A a {s) = Va~(K s (2a) + K s _ x {2a)) 
-ir(s)LT a {s) = -iB a (s) = Ml<s(2a) - K s _x{2a)) 

poo 

A a (s) - iB a (s) = £ a (s) = 2ja~ K s {2a) = y/E / e -«(*+i)^-i dt 



A a (s) 
B a {s) 



They verify the first order system, where fi(a) = a<p + (a) + a<p (a) = 2a: 



1 

-1 



da 



fx(a) 
fi(a) 



A a (s) 

B a {s) 



<*-h 



(176g) 
(176h) 

(176i) 



(176j) 



The pair 



A a (s) 

B a (s) 



is the unique solution of the first order system which is square-integrable with 



respect to eHog(a) at +oo. The total reflection against the exponential barriers at log(a) 



of the associated Schrodinger equations realizes 
matrices. 



r(i-5 
rfTj 



and 



r(i-s) 
r( 8 ) 



(1ft(s) = t:) as scattering 



From (163) we have A a {\) = \/ir e~ 2a . To normalize A a according to A a {\) = 1, we would have 
to make the replacement A a — > tt - ^ e 2a *4 a and B a — > y^TT e~ 2a B a and the expression of £ a in terms 
of the A'-Bcsscl function would be less simple. Let us also note that according to (116) wc must 



have g - l(2Q) 
nave Ka (2a) 



Regarding the isometric expansion, as given in theorem 15, we apply it to a function F(s) = 
/ °° k(x)x~ s dx such that -^xk(x) as a distribution on R is in L? . Using the L 2 -function ^A a (s), 
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which is the Mcllin transform of the function C a (x) = ^ A a (x) dx, and the Parseval identity we 
obtain a(u) = 2 f °°(— x-^k{x))C a {x) dx as an absolutely convergent integral. The square integrable 
function C a (x) is explicitely: 



C a {x) = ^l x>a (x)(j (2^/a(x-a)) + y Ji(2V 'a(x - a))) (177) 

And under the hypothesis made on x-^k(x) we obtain the existence of: 

f x d 

a{u) — lim y/ak(a) — 2Xk(X)C a (X) + yfa \ k(x)(l + —)J (2^/a(x - a)) dx (178) 

X~ >oo J a OX 

Let us observe that k(x) = ^ dy with l(y) g L 2 , so \k(x)\ 2 < f and Xk{X)C a {X) = 0(X~i). 
Hence: 



a(u) = Vafc(a) + V« / fc(x)(l + — ) J (2 ^fa(x - a)) (179) 

Comparing with equation (20a) we see that the f(y) defined there is related to a(u) , u = log(a) by 
the formula f(y) = a (log (a)), a = |, so \f(y)\' 2 dy = ^|a(log(a))| 2 2 da = |a(log(a))| 2 ±dlog(a). 
Similarly we obtain (3{u): 

d 



(3(u) = y/Ek(a) + y/a J k(x)(-l + — )J (2y/a(x - a)) dx (180) 

and comparing with (20b), g(y) = ^^/3(log(a)), \g(y)\ 2 dy = |/3(log(a))| 2 idlog(a). So according to 
theorem 15 we do have equation (20d): 

/>oo 

(\f(y)\ 2 + \g(y)\ 2 )dy= / \k(x)\ 2 dx (isi) 

Jo 

From 15 the assignment k — > (a, (3) extends to a unitary identification L 2 (5R(,s) = |; ■^■)~L 2 (M — > 
C 2 ; 4^), which has the property 7i(fc) — > (a, — /3). In order to complete the proof of the isometric 
expansion, it remains to check the equation (20c) which expresses k in terms of / and g. According 
to 15 we recover k(x) has the inverse Mcllin transform of J R (a(u) 2A a (s) + fi(u) 2(— iB a (s))) 4p. 
Expressing this in terms of f(y) and g(y), y = 2a, u = log(a), this means the identity of distributions, 
where we suppose for simplicity that f{y) and g(y) have compact support in (0, +oo) (as usual, this 
means having support away from as well as oo.): 



k{x) 



=r ( 2 vf /(y)2A * (a;)+2 \/i 3(2/)2H ^ (x)) ) % 

/ y (182) 

= 2 / (Vyf(2y) Ay(x) + ^g(2y) 2{-iB y {x))) 



Then imagining that we are integrating against a test function tp(x) and using Fubini we obtain: 

2 J°° ^f(2y)A v {x)^ 

poo / i \ (I 83 ) 

= J f(2y) U(x -y) + l x>y (Jo(2^y(x - y)) - J -H_J 1 (2^ - y))J J dy 
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(184) 




Proceeding similarly with g(y) one obtains for 2 J °° *Jyg(2y){— i B y (xj) — : 




J (y/y(2x - y)) + 



2x-y 



V 



J 1 (^y(2x-y)))f(y) dy 



(185) 



Combining (183) and (185) in the formula (182) for k(x) we obtain equation (20c). 

8 The reproducing kernel and differential equations for the 
extended spaces 

Let L a C i 2 (0, oo; dx) be the Hilbcrt space of square intcgrablc functions / which are constant in 
(0, a) and with their 7i-transforms again constant in (0, a). The distribution x-^ -^x f = -jj^x %-rzf 
vanishes in (0, a) and its Ti transform does too. So s(s — l)f(s) is an entire function with trivial 
zeros at — N. The Hilbcrt space of the functions s(s — l)T(s)/(s) satisfies the axioms of [2]; we prove 
everything according to the methods developed in the earlier chapters. Our goal is to determine the 
evaluators and reproducing kernel for the spaces L a . 

For / € L a , f(s) is a meromorphic function with at most a pole at s = 1 and also /(0) 
does not necessarily vanish. The Mellin-Plancherel transform f(x)x~ s dx = c(f)x~ s dx + 
f°° f(x)x~ s dx has polar part — jzi- Let us write (/, Y±) = — c(/) = Res(/(s),s = 1) = s(s — 
l)r(s)/(s)| s= i. This defines an clement Y x e L a . We define also ^ = r(l)Fi = Yi. Then (/,^f) = 
s(s - l)T(s)f{s)\ s =i. We also define y§ as the vector such that (/,^o) = s(s - l)T(s)f(s)\ s=0 = 
-/(0). One observes (f,H(Vi)) = {H{f),Y x ) = s(s - l)T(s)H(f)(s)\ s=1 = s(s - l)T(l - s)f(l - 
s)\ s= i = — /(0) = (f,y§) so yg = H(yi). To lighten the notation we sometimes write y± and y$ 
instead y x and y$ when no confusion can arise. 

We will also consider the vectors Xf € L a such that V/ G L a f(s) = (X* , /). 30 The orthogonal 
projection of X* to K a C L a is X s . Let us look more closely at this orthogonal projection. 
First let iV a be the (closed) vector space sum L 2 (0,a;dx) +HL 2 (0,a;dx). Inside N a we have the 
codimension two space M a defined as the sum of (lo<x<a)^ H L 2 (0,a;dx) and of its image under 
H. Finally, let R a be the orthogonal complement in N a of M a , which has dimension two. For a 
function / to belong to L a it is necessary and sufficient that its orthogonal projection to N a be 
perpendicular to the functions in L 2 (0, a; dx) which are perpendicular to lo<x<a, and the same for 
the 7i-transform, so this means exactly that its orthogonal projection to N a belongs to R a . So we 
have the orthogonal decomposition of L 2 (0, oo; dx) into the sum of the three spaces K a , R a and M a 
and L a = K a © R a . For / 6 L a to be in K a it is necessary and sufficient that c(f) = — (/, y? ) = 
and the same for c(H(f)), so this means that {3^1, 3^} is a basis of R a - The function y± belongs 

30 sometimes written X$ x . 
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to N a = L 2 (0, a; dx) + HL 2 (0, a; dx) and as such is uniquely written as u\ + Hvi. As yf £ L a we 
have constants a, /3 £ C such that: 



u\ + H a vi = -a 
H a iii +vi = —j3 



(186a) 
(186b) 



where we recall that P a is the restriction to (0, a) and H a = P a HP a , D a = H 2 . From what was said 
previously a = (y?,yf) and (3 = (H(y?),y?) = (y§,yf). We have thus: 



ui = (1 - D a ) 1 (-al 0<x<a + (3H a (l 0<x<a )) 

V! = (1 - D a )- 1 {+aH a {l 0<x<a ) - (31 0<x< a) 



(186c) 
(186d) 



Also the function y% may be obtained as the orthogonal projection to L a of — ^lo<a;<a- Indeed 
it follows from what has been seen above that for any element / £ L a , (f,y?) = — \ J Q a f(x) dx. As 



the function 



L 0<:r<a 



already belongs to N a , we have 



1 

— lo<a:<a 

a 



ui + Hvi + u 2 + Hv 2 



(187) 



where u 2 + Hv 2 belongs to M a , which means that u 2 £ L 2 (0 7 a;dx) verifies J Q ° u 2 (x) dx = and 
v 2 £ L 2 (0, a; dx) verifies v 2 (x) dx = 0. But there is unicity so we have exactly 



Ml 



U-2 — lo<x<a 

a 



V\ + v 2 = 



And we deduce: 



pa pa 

i u\{x)dx=—l I vi(x)dx = 

Jo Jo 

So a and (3 are determined as the solutions of the system: 

a (lo<x<a, (1 — D a ) 1 lo<x<a) — P(lo<x<a, (1 ~ D a ) 1 H a lo <x<a ) = l 

±0<x<a ±0<x<a ) = o 

We thus have: 
Proposition 22. Let p(a) and q(a) be defined as 

P(a)= [ a (l-D a )-\l 
Jo 

q(a)= [ (l-£> a )-\ff a (l 
Jo 

then: 



(188) 
(189) 

(190a) 
(190b) 



t<a )(x)dx 
o<x<a){x)dx 



p(a) -q{a) 
-q(a) p(a) 



{y ,yi) (yo,y a 





"1 


0" 







1 



The evaluators y± and y$ = 7i(3^i ) are given as u\ + Hvi and 7iu\ + v\ with: 

ui = -(y l5 yi)(i - D a )-\i 0<x<a ) + (y ,yi)(i - ^a)" 1 ^!^^) 
vx = -(y ,yi)(i - D a )-\i 0<x<a ) + (y Q ,y )(i - d^h^Uk*^) 



(191a) 
(191b) 

(191c) 

(191d) 
(191e) 
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We have introduced, for s ^ 0,1, X* as the evaluator f(s) for functions in L a . We shall write 
X* = T(s)X* and then y s = s(s — ■ This is compatible with our previous definitions of yf 
and ■ We note that the orthogonal projection of X* to K a is X s . So we may write 



x* =x s + \( s )y? + fi(s)ys 



(192) 



We shall also write y a (s,z) = / ° c ' y*(x)y%(x) dx = z(z - l)T{z)y}(z). One has H(y a s ) = y?_ s 
and y a (s, z) = y a (l — s, 1 — z) = y a (z, s). Taking the scalar products with y" and y§ in (192) we 



obtain 



1 



8(8-1) 



y a (l,s) = \(s)a + »(s)0 



1 



X(s) 
/i(s) 



8(8-1) 
1 



y a (i,i 



\(s)/3 + fx(s)a 



8(8-1] 
1 



( P y a (i,s)- q y a (i,i-s)) 
(-qy a (i,s)+ P y a (i,i-s)) = x(i-s) 



(193a) 
(193b) 

(193c) 
(193d) 



s(s - 1) 

Combining with (192), this gives: 

y s = s( s - i)x s + y a (i,s)( P y? - q ys) + y a (i, 1 - s)(- q y?+ P y$) (194) 

Let T a (s) = P (a)y a (l, s) - q(a)y a (l, 1 - a) (195) 

Proposition 23. The (analytic) reproducing kernel y a (s,z) of the extended space L a is given by 
each of the following expressions: 



s( s -i)z(z-i)x a (s,z)+ [y a (i,s) y a (i,i-s)] 

= s(s - l)z(z - l)X a (s, z) + [T a (s) T a (l - s)] 

= s(s - l)z(z - l)X a (s, z) + T a (s)y a (l, z) + T a (l - s)y a (l - z) 



p(a) 


-q(a) 




y a (i,z) 


-q(a) 


p(a) _ 




y a (i,i- 




a(a) 


0(a) 






Ta(z) 




0(a) 


a(a)_ 




T a (l - z) 





(196a) 

(196b) 
(196c) 



A very important observation, before turning to the determination of the quantities p(a) and 
q(a) shall now be made. Let L be the unitary operator: 

L(f)(x) = f(x)--[ X f(y)dy (197) 
x Jo 

It is the operator of multiplication by at the level of right Mellin transforms. Obviously it 
converts functions constant on (0,a) into functions vanishing on (0, a). Let us now consider the 
operator 

H° = LH L- 1 = L 2 H = HL- 2 (198) 
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It is a unitary, self-adjoint, self-reciprocal, scale reversing operator whose kernel is easily computed 
to be 

It has L(e~ x ) = (l + ^)e~ x — - as self-reciprocal function; the Mcllin transform is ^-^r(l — s) which, 
multiplied by s(s — 1) gives (1 — s) 2 T(l — s) which is the Mellin transform of a more convenient 
invariant function for 7i°, the function x(x — l)e _x . This function is the analog for Ti of eT x for 
Ti. Let us now consider the space L(L a ). It consists of the square integrable functions vanishing 
identically on (0, a) and having Ti® transforms also identically zero on (0, a). But then the entire 
theory applies to Ti exactly as it did for Ti, up to some minor details in the proofs where the function 
Jo was really used like in Lemma 9 or proposition 11. We have for Ti versions of all quantities 
previously considered for Ti. To check that the proof of 11 may be adapted, we need to look at 
K{x)= / *fc°(i)dt = V^i(V^)+2(Jo(2V^)-l) + 2 ^kMH) rfjl=0 ^ +oo(l i). So we may 
employ lemma 8 as was done for Ti. Proposition 11 and Theorem 12 thus hold. We must be careful 
that the operator L~ l is always involved when comparing functions or distributions related to 
with those related to Ti. For example, one has Xf = -^L~ l X^ and y s = s(s - l)X* = L~ X X^. 
The two types of Gamma completed Mellin transforms differ: for Ti we consider T{s)f(s) while for 
Ti° we consider s 2 T(s)g(s). Indeed this is quite the coherent thing to do in order that: 

s 2 T(s)g(s) = S ( S -l)r( S )f( S ) (200) 

for g = L(f). The bare Mellin transforms of elements of spaces K° are not always entire in the 
complex plane: they may have a pole at s = 0. After multiplying by s 2 F(s) which is the left 
Mellin transform of the self-invariant function x(x — l)e~ x , as F(s) is the left Mellin transform 
of e~ x , we do obtain entire functions, whose trivial zeros are at — 1, —2, ...(0 is not a trivial 
zero anymore.) From equation (200) we see that the (analytic) reproducing kernel X°(s,z) exactly 
coincides with the function y a (s 7 z) whose initial computation has been given in Proposition 23. 
Also the Schrodinger equations will realize ± (^ s -) 2 ^"fey^ as scattering matrices, and there will be 
an isometric expansion generalizing the de Branges-Rovnyak expansion to the spaces L a . We will 
determine exactly the functions A^(s), B°(s), £ai s ) an d especially the function /i°(a). It will be 
seen that this is a more complicated function than the simple-minded fi(a) = 2a. . . 

The key now is to obtain the functions p(a) and q(a) defined in Proposition 22, and the function 
y a (l,s). It turns out that their computation also involves the quantities (we recall that Jq{x) = 
J (2V^)): 



r(a)=l+ ((l-Da^Ha-J^dx (201a) 
Jo 

s(a)= [ ((1-Da)- 1 ■ JS)(x)dx (201b) 
Jo 

In order to compute r, s, p, q we shall need the already defined functions 4>t (= (1 + Ha)^ 1 JS on 



47 



(0, a)), 4> a , (= (1 — H a ) 1 Jq) as well as the entire functions ^+ and V> a verifying: 

tJj+ +HP a i'+ = 1 (202a) 
V" - HPair = 1 (202b) 

We have r(o) = 1 + ^ / a (— 0^ (x) +^>~(x)) dx, and we know explicitely <jr£. But, we shall proceed in 
a more general manner. First we recall the differential equations (12fa), (121b) which are verified 
by 4% (where S x = x-§^ + §): 

a lk^ = +6x ^ ~ (A<(a) + \^ (2 ° 3a) 
= +Mj + (A*(o) - (203b) 

We compute ar'(a) = a ' " (a) 2 +fc(a) +i/ a (a:^ + l)(^(s) - ^-(x))+/i(a)(<^.+ (x) + ^-(a;)) and 
simplifying this gives exactly ar'(a) = £t(a)| fo{4>£ ( x ) + ( fia( x ))dx = s(a). Similarly starting with 
s ( a ) = kfoi&tW+fc (x))dx we obtain as' (a) = \n{a) + \ f£ (x (x)+(j)- (x))+n (a) (-(/)+ (x) + 
4>a( x ))) dx which gives fi(a)r(a) — s(a). So the quantities r and s verify the system: 

ar'(a) = (j,(a)s(a) (204a) 
(as)' {a) = /x(a)r(a) (204b) 

Either solving the system taking into account the behavior as a — > or using the explicit formulas 
for (/>£; we obtain in this specific instance of the study of Ti, for which fJ,(a) = 2a, that r(a) = Iq(2o) 
and s(a) = I\(2a). 

From (202a) and (202b) we obtain two types of differential equations, either involving or 
a&. From ^+(x) + /„" J (2^xy)i>+ (x) dx = 1 , we obtain (1 + HP a )a^+(x) = -a^+(o) J a . We 
do similarly with tjj~ and deduce: 

a ^t( x ) = -a*l>i(.a.)<t>i(z) (205a) 
a-^ipa (x) = +ai'~ (a)<t>- (x) (205b) 

Regarding the differential equations with x-^, which we shall actually not use, the computation is 
done using only the fact that the kernel is a function of xy so xJ^Jo{2y/xy) = y-g-Jo(2- s /xy). We 
only state the result: 

{X ~&x~ + \ ] ^ {X) = \^ {X) ~ ^M^) (206a) 

+ \ ] ^ {x) = ^aW+a^ia^tix) (206b) 

Let us now turn to the quantities p(a) and q(a). We have p(a) = J a (l — Da)^ 1 {lo<x<a)(x) dx = 
U a (^(x)+ip-(x))dx. Sop'(a) = ^+(a)+i>-{a))-±i>+(a)fZcf>+(x)dx+^- (a) ft <(>- (x) dx. 
Reorganizing this gives: 

p , {a) = V>+(o)+^(o) (1 + r -0+(x)+0-(x) ^ + -^+(a)+^-(a) r +0+(.t)+0-(x) ^ 
2 Jo 2 2 Jo 2 

(207) 
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We remark that from the integral equations defining ip^ we have ipa( a ) = 1 — Jo Jo(2\ // a^)?P£ (%) dx = 
1 - /„" cf>+(x) dx and 1>~ (a) = 1 + /„" J (2vM^a (*) = 1 + Jo" <fc (*0 So ^" +(a) + fc(a) = r(o) 
and ^» ( - a ^ +,;: '° ^ = s(a). Hence the quantity p(a) verifies: 

p 1 '(a) = r(af + s{a) 2 (208) 

With exactly the same method one obtains: 

q'(a) = 2r(a)s(a) (209) 

Let us observe that q(a) = \ f£((l - Ha)- 1 - (1 + H a )- l )(l) dx = 0(a 2 ) and p{a) = \ / Q a ((l + 
H a )~ l + (1 - Ha)- 1 )^) dx — a ^o a. So (p ± q) ~ a ^ a. Also r(a) ~ a ^o 1 and s(a) — a ^o a- The 
equation for p(a) can be integrated: 

p(a) = a(r(a) 2 - s(a) 2 ) (210) 

Indeed this has the correct derivative. Regarding q(a) the situation is different, one has q' = 2rs = 
^rr' so in the special case considered here, and only in that case we have q(a) — \(r(a) 2 — 1). 
Summing up: 

Proposition 24. The quantities r(a), s(a), p(a) and q(a) verify the differential equations ar'(a) = 
/j,(a)s(a), as'(a) + s(a) = n(a)r(a), p' '(a) = r(a) 2 + s(a) 2 , q'{a) = 2r(a)s(a), p(a) = a(r(a) 2 — s(a) 2 ). 
In the special case of the H transform one has: 

r(a) = I (2a) (211a) 

s(o) = 7i(2a) (211b) 

p(a) = a(I${2a) - if (2a)) (211c) 

q(a) = ±(I 2 (2a)-l) (211d) 

We now need to determine y a (l,s) — s(s — l)T(s)Yi(s). There holds Y\ = ii\ + Hv\ — 
-alo<i«. + lxyaHvi- So Y 1 (s) = a 2 ^- + (Hvi)(x)x- S dx. Then $™ (Hvi)(x)x~ s dx = 
J °° vi(x)g s (x) dx = Jq Vi(x)g s (x) dx, where the function g s from (73) has been used. Recalling 
from (76a), (76b) the analytic functions u s , equal to —(1 — D a )- 1 H a (g s ) on (0,a), and v s , equal to 
(1 — Da)- 1 P a (g s ) on (0, a), and using (186d) and self-adjointness we obtain 

/•a pa 

V\(x)g s (x) dx = ~a I u s (x)dx — [3 / v s (x)dx (212) 
Jo Jo 

Let us now recall that we computed ((83a)) [x-g^ + s)u s and found it to be on the interval (0,a) 
given as — av s (a)(l — D ) _1 (Jo) — a(a~ s + u s (a))(l — D a )~~ 1 H a {J 1 )- Integrating and also using 
equations (99) and (102) we obtain 

Va~Ea{s) - a 1 ' 8 + (s - 1) / u s (x) dx = -^iH(E^)(s)s(a) - Va^(*)(r(o) - 1) (213) 



n 
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u s (x) dx = \fa 



a—* - E a (s)r(a) -H(E a )(s)s(a) 



s - 1 



(214) 



We have similarly ((83b)) ( x £+l- s ) Vs = -^EE^l-D^i J^-y/EHiE^a^l-D^H^ Jg) 
so integration gives av s (a) — s v s (x) dx = — ^/aE a (s)s(a) — yfaTt(E a )(s)(r(a) — 1) hence 



v s (x) dx = \fa 



E a (s)s(a)+H(E a )( S )r(a) 



(215) 



Combining (214), (215) with (212), and using y a (l, s) = s(s - l)T(s)Yf(s): 

'a(a)r(a) f3(a)s(a)\ /-^j7n~\t \{ a ( a ) s ( a ) P( a ) r ( a ) 



Y 



(s) = s/aE a (s)[ - 

V s — 1 



) (216a) 



y a (l, s) = ^E (sa(a)(£ a (s)r(a) + £ a (l - s)s(a)) + (1 - s)p(a)(£ a (s)s(a) + £ a {\ ~ s)r(a))) 

(216b) 



Proposition 25. The functions 3a(l, s) and y a (l, 1 — s) verify 



a(a) /3(a) 
[3(a) a(a) 



s(£ a (s)r(a) + £ a {l - s)s(a)) 
(l-s)(£„(«)s(a) + £ Q (l-s)r(a)) 



(217) 



Comparing with equation (195) we get: T a (s) = y /r as(£ a (s)r(a) + £ a (l — s)s(a)). So: 
Theorem 26. T/ie analytic reproducing kernel y a {s^z) associated with the extended spaces L a is: 



y a {s, z) = s(s ~ l)z(z - l)X a (s, z) + [T a (s) T a (l - s)] 



a(a) /3(a) 
/3(a) a(a) 



Ta(z) 

T a (l - z) 



(218a) 



X a (s, z) 

Ta(s) 

a(a) 



£ a (s)£ a (z)-£ a (l-s)£ a (l-z) 



£ a (s) = 2y/EK a (2a) 



(218b) 



s + z-1 

T a (s)=^s(£ a (s)r(a)+£ a (l-s)s(a)) r(o)=/ (2o) s(a) = h(2a) (218c) 

p(a) 



(3(a) 



p(a) 2 - q(a) 2 

g(o) 
p(a) 2 - q(a) 2 



p(a) = a(I 2 (2a) - It(2a)) 
g(a) = l(/ 2 (2a)-l) 



(218d) 
(218e) 



We proceed now to the determination of A%, B° and = A^(s) — iB^(s). The function A^(s) 
is even under s — > 1 — s and S*(s) is odd. We must have: 

zy a (l, z) = 2(-iBZ(l))Al(z) + 2At(l)(-iBt(z)) (219) 

On the other hand from (195) we have y a {l, z) = aT a (z) + (3T a (l — z). Let us write 

zT a (z) = y/a(z(z - l)r(a)£ a (z) + zr(a)£ a {z) 

+ z(z - l)s(a)£ a (l - z) + zs(a)£ a {\ - z)) (220) 

zT a (l -z) = VE{-z{z - l)s(a)£ a (z) - z(z - l)r(a)£ a (l - z)) (221) 
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zy a (l, z) = V^(z(z - l)((or - (3s)£ a {z) + (as - (ir)£ a {l - z)) 

+ zu(r £ a {z) + s£ a {\ - z))} (222) 

Extracting the even part {zy a (l, z)) + and the odd part (zy a (l, z))~: 

(zy a (l, z))+ = V^(z(z - l)(a - (3)(r + s)A a + (z - ^)a{r - s){-iB a ) + ^a{r + s)A a ) (223) 

(zy a (l, z)y = V^(z(z - l)(a + /?)(r - «)(-»B«) + (z - ^)a(r + s)A a + ^a(r - *)(-*B„)) (224) 

We have («y a (l,«))+ = 2(-i6*(l))^(z) and (aO^l,*))- = Z4*(1)(-*B£(;»)). Let us define 
K(a) = (2(-iB*(l)))" 1 and L(a) = (2.A*(1)) _1 . We know that: 

= £r£ = i ( 225 ) 

So it must be that 

K(a)(a - 0)(r + s) = L(a)(a + (3){r - s) (226) 

Also, taking z = 1 in (223) we have = 2y/a\a (r£ a (l) + s£ a (0)) = aT a (l). But referring to 
(195) one has T (l) = pa - q(3 = 1. So: 

A»L(a) = -L (227) 
a(a) 

Then: 

K( ,2 = 1 (« + /3)(r-.s) a 2 ~/3 2 r 2 - s 2 lg 1 

W a(a)(a-0)(r + s) a(a) (a - /3) 2 (r + s) 2 p a (a - /3) 2 (r + s) 2 1 j 

{a- /3)(r + s)K(a) = a _ 3 (229) 

We conclude: 

At = z(z l)A a + (z- \) a[r R 7, S \ A -iB a ) + (230a) 

2 (a — p)(r + s) l(a — p) 

-iBl = z(z - l)(-iB ) + (z - 1) " ( ^+ 8) A + — ^r(-iB a ) (230b) 

2 (a + pj(r — sj 2(a + p) 

Let us now observe that — x« = -?- and further: 

a r-s p a(r-s) 2 p' - d 

-a— log(p-g) (231a) 



a — (3 r + s p — q p p — q da 

a r + s p a(r + s) 2 p' + q' d 



a- -=a— \og{p + q) (231b) 



a + (3 r — s p + q p p + q da 



A° a (z) = (z(z - 1) + lJL-)Aa(z) + a-^ Iog(p - g)(z - (232a) 
-iB»(z) = (z(z - 1) + I-^-)(-iB a (^)) + a A i og (p + 9)(z _ 1)^) (232b) 



Combining we get finally: 
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Theorem 27. The E function associated with the entire functions s(s — l)T(s)/(s), f £ L a is: 

Id..,.. 1. 1 



E*(z) = (z{z - 1) + -a— log(p(a) 2 - q(a) 2 )(z - -) + -p(o)a(o)) 

/l d , p(a) + q(a) 1 1 , ... . 



(233) 



where p(a) = a(/ 2 (2a) - if (2a)), q(a) = ±(J 2 (2a) - 1), a(a) = p(a ff (a) , , /3(a) = p(a ff (n) * , and 
= 2V^A" z (2a). 

We shall now obtain by two methods the function [J? (a). First, we compute £a(h) = (— \ + 
\p{a + p))£ a {\) = \*%£ a (\) and invoke a&S>{\) = -fi«(a)£°(±). We thus have: 

Theorem 28. The mu function for the chain of spaces L a , < a < 00 is 

l f{a) = n{a)+a±\og 1 ^_ (234) 
da p + q 



2a + a±lo g ? a - 1 \i? a ]- 2a f? a ] + 1 (235) 
da 8 (2a + l)/ 2 (2a) - 2al\ (2a) - 1 1 ; 

2a- 2 + o(l) (a ^+00) (236) 



The asymptotic behavior is a corollary to lim a ^oo 2a pjpr = —2 which itself follows from 
p 2 — g 2 ~ |;/o(2a) 4 Y^T and (-)' ~ + which are easily deduced from the asymptotic expansion 
I (x) = ^==(1 + h + S? + •■ ■) (I 33 ])- Of course the o(l) is in fact an 0{a- r ). 

The second method to obtain //*(a) relics on £a ^}^ > ^a-^+oo ((116)). We have: 

£ » _ 1+ MM^!M +0( i } (237a) 



a 2 £ a (a) a a 2 

£lil - a) - X-JcaA^Eli) 2 (237b) 



er 2 £ a (l — cr) 2 da p — q At(a) 

so ^ff = 1 ~ JIJ7I) a lk 1°S an< i (234) is confirmed. We can use this method to gather more 
information. From (115a) we have, as 5R(s) — > +00: 



+ ^^^ + 0(l)) 
- «*-(! + ^ +(a) 2 7^ (a) + 0(1)) (238b) 



Let us be careful that £ a (s) = T(s)E a (s) while £*(s) = s 2 r(s)^(s). We obtain: 

a0 o+ (a) - a0°~(a) = a0+(a) - a0 _ (a) + a-f- log(p 2 - g 2 ) - 2 (239a) 

da 

a<£° + (a) = a0+(a) + a 4" log — - (239b) 
aa a 

a^- (a) = a0- (a) - a4- log (239c) 
aa a 
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We recall that (p± q) ^ a _»o a - We integrate (239b) and (239c) using (130a), (130b) and this gives 
det(l + H°) = ^ det(l + H a ) and det(l - H%) = £±2 det (i _ Ha ). 

det(l + HT) = ^—^ det(l + H a ) = dct(l + H a ) - / (r - s) 2 da (240a) 
a a J 

det(l - HT) = det(l - H a ) = det(l - H a ) - f (r + s) 2 da (240b) 
a a J 

Theorem 29. Lei 7i° = L7YL -1 be the self-reciprocal operator on L 2 (0 , oo; dx) with kernel: 



, M2^Sy) l-J (2^y) ^ n 2 x n y n 

J {2y/xy)-2 f=— + * = > -1 ? — — - 

y/xy xy (n + l)r 

and let H% be the restriction to L 2 (0,a;dx). Then: 



(241a) 



det(l + HT) = e +a -i a " - [" \l (2a) - h(2a)) 2 da = e +a '^ (l 2 (2a) - I 2 (2a) - - ) 
a Jo 2a ' 

(241b) 

2a ' 
(241c) 



det(l - HT) = e - Q -^ 2 - f (I Q (2a) + h(2a)) 2 da = e~^ a ' (/ 2 (2a) - L?(2a) 
a Jo 



From theorem 26 ||^i || 2 = ±\\Xi \\ 2 + 2(a + (3)T a {\) 2 and T a (±) = \y/a(r + s)E a {\), Also, 
(a + /?)(r + ,s) 2 = -i-(r + S ) 2 . Furthermore ||yi|| 2 = ^(i) 2 !) A? || 2 = ^||A?|| 2 and j|Ai|| 2 = 
7r||A.ij| 2 . And also from (163) £ a (§) = dct(i+g") anc ^ f rom theorem 17 one has ||A?|| 2 = 
2 J a °° (det fe^) 2 f and the analog holds for XT. Let us observe that XT = —LX* so ||A? a || = 



A? x ||. So 

2 

nvaxiis „ f°° ( ', l-H?\ 2 db n [°° ( , l-H b \ 2 db n p' + q'f, 1 -# a x 2 

11*4*11 = 2 / det-— | -=2/ det--/ - + 8a 1 det ■ 



i + #?y & 7 Q V l + H b ) b p-q \ l + H, 



(242) 



Theorem 30. Let L a 6e t/ie Hilbert space of square integrable functions on f G L 2 (0, oo; da;) suc/i 
£/iai 6oi/i / and 7i(J) = J °° Jo(2yfxy)f(y) dy are constant on (0, a). TTien i/ie squared norm of the 
linear form f i— ► J °° dx is given by either one of the following two expressions: 

(2b + l)r m -2bl 2 { 2b)-l^e^ db (243&) 



(26 - l)/ 2 (26) - 2bl 2 {2b) + 1J b 

= 2 r db + 2 _ ^(±^A ^ (^b) 
y a 6 (2a-l)Ig(2a)-2aZ?(2a) + l v ; 

The squared norm of the restriction of the linear form to the subspace K a of functions vanishing on 

(0, a) and with H(f) also vanishing on (0, a) is 

One may express the wish to verify explicitely from equations (230a) and (230b), or in the 
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equivalent form 

Al(z) = ((z - \) 2 + i£±|) A a (z) + {aj- a Iog(p - g))(* - ^(-iS.C*)) (244a) 
-i^(z) = ((z - i) 2 + ||^|) (-AB.(*)) + (o^ log(p + - i)A(z) (244b) 



the differential system: 



a^^(z) = - M »^(z) - (z \){-iK{*)) (245a) 
a^(-^( Z )) = +f(a)(-iBl(z)) -(z- l -)Al{z) (245b) 
and also to verify explicitely the reproducing kernel formula 

y ais , z) = mm^^m^i (246) 

The interested reader will see that the algebra has a tendency to become slightly involved if one 
does not benefit from the following preliminary observations: using p' = r 2 + s 2 , q 1 = 2rs, ar 1 = fir, 
as' = fir — s, p = a(r 2 — s 2 ) one first establishes aq" + q' = 2fip', ap" + p' — 2 = 2fiq'. Using this 



one checks easily: 



Also the identity 



P + 1 J V P + 1 J a z p + q a p + q 

pwv + fpwy = i_p _ (247b) 

p-q J \p-q J a-p-q a p-q 

p 2 p' + q' p' - j_ , nA ~ \ 

= a a = pa (247c) 



p 2 — q 2 P + q p — q 
is useful. The verifications may then be done 



9 Hyper functions in the study of the TL transform 

In this final section we return to the equation (31): 

W)(i*) = ^7(*'^), (248) 
Let us recall that / £ L 2 (0,oo;dx) and ip : L 2 (0,oo;dx) — > L 2 (0,oo;dx) is the isomctry which 
corresponds to F(w) h-> F(w 2 ), where F(w) = X^o ™ - " 7 "' f( x ) = 2^=0 c n P n (x)e~ x , P n {x) = 
Ln (2x). Let g = tp(f). Using A = it, in the L 2 sense: 

-i />oo \ I • \ 1_ 

^sL^-^^^ (249) 

It is natural to consider separately A > and A < 0. So let us define: 

i r° A+i~A-i 



G+{X) = 2^]_ -2T f{ 2 )e dX (25 ° a) 
G-(x) = ±±Lf(±zA)e-^d\ (250b) 
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We observe that G+ is in the Hardy space of > and G_ is in the Hardy space of 3(x) < 0. 
Their boundary values must coincide on (—00, 0) as g £ £ 2 (0, +00; dx). So we have a single analytic 
function G(z) on C \ [0, +00) with G = G+ for S(x) > and G = G_ for S(x) < 0. Then 
g = ip(f) = G+ — G- is computed as 

g(x) = G(x + iO) - G(x - iO) (251) 

In other words g is most naturally seen as a hyperfunction [23] , as a difference of boundary values 
of analytic functions. We shall now compute it cxplicitcly, and also we will show later that this 
observation extends to the distributions A a (x), —iB a (x), E a (x) which are associated with the study 
of the Ti transform. The point of course is that the corresponding functions G will for them have a 
simple natural expression. 



We have, for S3(z) > 0: 



-1 poo \ ; \ I 1 

G( Z) = H ^R-^^dX (252a) 



G ( z ) = I e^- x+ ^f(x) dx ) e+* Xz dX (252b) 



2tt J 2A VJo 



Let fi = |(A — j), X = fi + \Jl + (X 2 , ^ = I+i with, for < A < 00, —00 < // < 00. 

g{z) = h r att or da; ) e+az ^ (252c) 

For /i — * +00, A = 2/i+^- + . . . , and for /i — > —00, A = — ^7 + - ■ ■ , and ~ ^- as — > —00. So far 
the inner integral is in the L 2 sense. We shall now suppose that / and /' are in L 1 (so lim-z^oo fix) — 
0) and write e~^f(x) dx = e~^-*e*f(x) dx = + ^ JjJ» e -*^(/(i) + /'(x)) dx. 

1 / r 00 x e +iXz r 00 x e +iXz ( r°° ■ \ 

2?r V A + z t/i + 1 J_ 00 X + 1 ifi + 1 \J J 

(252d) 

In this manner, with / £ L 1 , f £ L 1 , ^( z ) > 0; we have an absolutely convergent double integral. 

G(z) = — (f(p) d^+ / ^ (/(x) + /'(x))rfx 

(252e) 

Observing £ ^ £±£ d M = £ / °° 4^-L_ e +'^ dA = ^ / °° ^e+ iA * dA, we then suppose 
5i(z) < 0, 3(z) > (or 3(z) > 0) so that we may rotate the contour to A = —it, < t < 00. This 
procedure gives thus: 

1 f°° X e +lXz , 1 f°° e tz 



d » = 7T- T—dt (252f) 



27r J.qc A + i ifj, + 1 2iri J 1 + t 

Also: 

■d f i = — J—e-tW-V*** d\ (252g) 



2 71 " J -00 A + z i/i + 1 27ri 
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We rotate the contour to A e z[0, — oo), which is licit as x > and, for SR(z) < 0, x > 0, we obtain: 

/ dt (252h) 

2mJ 1 + t y ' 

Going back this allows to write (252e), for $t(z) < 0, $s(z) > as: 

1 / r°° e zt r°° ( r°° e zt -§ (*+t) \ \ 

G < z »=^(«°'l — t it+ L U -T+i- < 252i > 

and finally, after integrating by parts: 

1 poo / POO 11 \ 

G{z)= 2rn] U 2 {1 + l )eZt ~ iV{t+i)dt ) f{y)dV (252j) 

This last expression (still temporarily under the hypothesis /, /' £ L 1 ) is certainly a priori absolutely 
convergent for $t(z) < and gives G(z) in this half-plane. 

We are led to the study of: 

i pOO 1 i 

a(z,y) = — -(1 + (253a) 
Z7TI j 2 i 

We still temporarily assume 5ft(z) < 0. We even suppose z < and make a change of variable: 

•<•■»> = 55 5 i"" * + 5 jf (2530 

o(2,! ' ) = 2S ylj^Tz K 'W«(y - 2 ")) + folVB - 22))) (253d) 
For any z £ C \ [0, +oo) and any y > the integrals in (253c) converge absolutely and define an 
analytic function of z. Furthermore the K Bcsscl functions decrease exponentially as y — > +oo 
in (253d). For fixed z, a{z,y) is certainly a square-intcgrable function of y (also at the origin), 
locally uniformly in z so the equation (252j) defines G as an analytic function on the entire domain 
C\ [0, +oo). Then by an approximation argument (252j) applies to any / <G L 2 (0, oo; dx) and any 
z e C\ [0,+oo). 

We now study the boundary values a(x+i0, y), a(x — i0, y), x, y > 0. We could use the expression 
of the K Bessel functions in terms of the Hankel functions and H^ 2 > , go to the boundary, and 
then recover the Bessel functions Jo and J\. But we shall proceed in a more direct manner. Let us 
first examine 

i i /.oo -1 

*?l J % , ' , r oo , f , f (254a) 



e 2 



2ni 2 J Q u 2ni J 1 ^Jt 2 - 1 

1 /"OO i />OC i -i 

_L / e -V^=^)* r i dt + _L / Wv(i>-a»)t ( _!: (254b) 

2m Ji 2ixi J 1 ^Jt 2 - 1 t 
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j -^y(y-2z)_ r°° e -Vvlv-**)*t- 2 dt 1 f 00 /T^T, 1 1 

= -= H + e-vub-^H . --)dt (254c) 

We now look at the (distributional) boundary values z — > x with z = a; + ie, e — > + or z = x — ie 
and e — > + . We shall take x > 0. Here the singularities at y = 2x and at y — arc intcgrablc and 
we need only take the limit in the naive sense. We distinguish y > 2x from < y < 2x. In the 
former case, nothing happens: 

d(x + iO, y) = d(x -i0,y) = — e -y/v(v-i*)t_^^ (25 5a) 

27TI J 1 V< 2 - 1 

In the latter case: 



1 +iy/y(2x-y) _ f°° +iy/y(2x-y)t f -2 J f i /-oo i -, 



d(x + »0,i/) = - p + — / e+'v 8 ^-^'; ; --)dt 

V ; 2tt ^ y (2x - y) 2m h V y/W^l * 

(255b) 

1 p -i\/y{2x-y) _ f 00 P -»V yi 2x ~y) * +-2 ^+ 1 /-oo -, 1 

V 2tt ^(2! - y) 27ri A V Vt^T t ; 

(255c) 

So + zO, y) — d(x — iO, y) is supported in (0, 2x) and has values there 



l cos^/y 2x-y cos(y/y(2x - y) t) t' 2 g 1 r 00 . _ 1 1 

j (n \ + - / sm(V2/(2x - y t){ - - (it 

(255d) 

We used this method to have a clear control not only of the pointwise behavior but also of the limit 
as a distribution. There is no necessity now to keep working with absolutely convergent integrals 
and we have the simple result, using the very classical Mehler formula: 31 



d(x + iO, y) - d(x - 20, y) = l < y <2x(y)- [ — — V ,^ X -- — dt = \ ±o <y< 2x{y)Jo{\f y{2x - y)) 

7T Ji y/t 2 - 1 2 

(256) 

Let us now consider the behavior of 

e ( Z) y) = -Li f e zt "^(*+T) dt = -L./— IL- i /°° e -f v^^X^) rfM (257) 
v yy 2™ 2 7 2™ V y - 2z 2 J y ' 

We make the simple observation that e(z,y) — J^d(z,y). So we shall have (as is confirmed by a 
more detailed examination): 

d 1 



e(x + iO, y) - e(x - iO, y) = — - l 0<y<2x {y)J (y / y{2x - y)) 

! j-g- j (258) 

S2x{y) - - lo< y <2x(y)J 2x _ J\W y(2x - y)) 



Combining all those elements we obtain that the function k(x) = ip(f)(x) is given as: 

k(x) = f(2x) + \ f Jo(\/v(2x-y))f(v) dy - \ I J 7 -^.hWy{2x-y))f(y) dy (259) 
2 Jo 2 J V 2x - y 



31 we are mainly interested in the boundary value as a distribution and we skip the discussion of the pointwise 
behavior at the borders y = and y = 2x. 
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Some pointwise regularity of / at x is necessary to fully justify the formula; in order to check if 
continuity of / at 2x is enough we can not avoid examining e(z, y) more closely as z — > x. 



(260) 



y -2z2m J 1 ^Jt 2 



2m y-2z y y -2z2TTiJ 1 ^/W^\ 



The integral term on the right causes no problem at all. And writing f 



y~2z ' y-2z 



- — V y~i z > ~ 1 ; again the term on the right has no problem, so there only remains , and of 
course, this is very well-known, the difference between +i0 and — iO gives the Poisson kernel, so for 
non-tangential convergence, continuity of / at 2x is enough. Of course this discussion was quite 
superfluous if we wanted to understand k as an L 2 function, here we have the information that non 
tangential boundary value of G(x + iO) — G(x — iO) does give pointwise the formula (259) if / is 
continuous at y = 2x. We can also rewrite (259) as: 

k{x) = (1 + Tx ] \ fj Jo ^ y{2x - y))f{y) dy (261) 

This is exactly one half of equation (20c), where k was obtained from (f,g) as ip(f) + w ' "0(ff)- 
Let us observe that w = verifies, as an operator, (^ + 1) • w = w ■ (^ + 1) = ^ — 1. 
So the isometry corresponding to g(w) i— > wG{w 2 ), which is the composite w ■ ip, sends g to 
(~ ^ 1x^1 Iq X ^o(^/y(^x — y))f{y) dy. This is indeed the second half of equation (20c). 

The formulas (20a) and (20b) may be established in an exactly analogous manner (taking k with 
compact support to simplify the discussion). But this would be a repetition of the arguments we 
just went through, so rather I will conclude the paper with a method allowing to go directly from 
<A a (s), ~iB a {s), £ a (s) to the distributions A a (x), —iB a (x), E a (x), and this will show that they arc 
in a natural manner (differences of) boundary values of an analytic function. 

From the expression £ a (s) = T(s)E a (s) = 2y/aK s (2a) = yfa J °° e~ a ( t+ ^H s ~ 1 dt, we shall recover 
E a (s) as a right Mellin transform with the help of the Hankel formula r(s)" 1 = j c e v v~ s dv, where C 
is a contour coming from — oo along the lower edge of the cut along (— oo, 0] turning counterclockwise 
around the origin and going back to — oo along or slightly above the upper edge of the cut. Let us 
write the Hankel formula as 

W) = ^il/ Vv ~ Sdv {t> 0) (262) 

So we have: 

Ea{s) = J o e tv v~ s dv^j e - Q ( t+ T) dt (263) 

Let us suppose 3?(s) > 1. Then the contour C can be deformed into the contour C e coming from 
— ioo to —ie, then turning counterclockwise from e~ l ^e to e l 2"e, then going to +ioo. Also we impose 
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< e < a. The integrals may then be permuted: 



E a {s) = ^E^-.j (f e tv e- ai - t+ ^ dt)v~ s dv (264) 



i) 



and using e(z, y) from (257) this gives: 

3?(s) > 1 => Ea(s)=Va I 2e(v,2a)v- s dv (265) 

We have previously studied e(z,y), which is also expressed as in (260). We see on this basis and 
simple estimates that we may deform C e into a contour C a ^ going from +oo to a + rj along the lower 
border, turning clockwise around a from a + r\ — iO to a + rj + iO, then going from a + rj to +00 on 
the upper border (rj <C 1). We will have in particular from (260) a term t^U J^^q 2 a-2v ^ v which 
is a _s . The final result is obtained: 

This identifies E a {s) as the right Mellin transform of the distribution 

E a (x) = y/a (s a {x) + l x>a (x)^-J (2y/a{x-a))^ = VE^- (lx>a{x)Jo(2y/a(x-a)j) (267) 

This proof reveals that the distribution E a (x) is expressed in a natural manner as the difference of 
boundary values y/a{2e{x + iO, 2a) — 2e(x — iO, 2a)). with 



1 r 1 / 

v^2e(z,2a) = — / e"*- ^ dfc = yfd~ —2 J K x {2Ja(a - z)) (268) 

2ni J a 2ni y a — z 

The formulas (176e) and (176f) are recovered in the same manner. 

Theorem 31. The distribution A a (x) = ^§{1 +H)(4>+l 0<x<oo ), (f>+(x) + f£ J (2^/xy)4>+ (y) dy = 
Jq(2^/ox) , is the difference of boundary values y / a(a(x + iO, 2a) — a(x — iO, 2a)), with: 



-1 l-OO 1 -1 

yfca{z,2a)=yfc — J -{l + - t )e zt - ai - t+1 -Ut 



'a 

2-ki 



(269) 



TTie distribution —iB a (x) = ^(-1 + W)(</> a lo<x<oo), </> a fa) ~ Jo" Jo(2y/x~y)(f> a (y) dy = J (2y/ax), 
is the difference of boundary values i/a(—ib(x + iO, 2a) — {—ib{x — iO, 2a))), with: 



1 l-OO 1 -1 

VE(-ib(z, 2a)) =V^^- i J o g (1 - i )e^- a ( t+ *) di 



2ttz 



~t ^^Xi(2V«(«-z)) - A' (2v/a(a - i))J 



(270) 
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10 Appendix: a remark on the resolvent of the Dirichlet 
kernel 



In this paper we have studied a special transform on the positive half-line with a kernel of a multi- 
plicative type k(xy), following the method summarized in [5, 6]. We have associated to the kernel 
the investigation of its Frcdholm determinants on finite intervals (0, a), and have related them with 
first and second order differential equations leading to problems of spectral and scattering theory. 
There is a vast literature on kernels of the additive type k(x — y), and on the related Frcdholm 
determinants on finite intervals. The Dirichlet kernel on L 2 (—s, s; dx): 

K s ( X ,y) = ^Liy) (271 ) 
n[x - y) 

has been the subject of many works (only a few references will be mentioned here.) The Fredholm 
determinant det(l — K s ), as a function of s (or more generally as a function of the endpoints of 
finitely many intervals), has many properties, and is related to the study of random matrices [22]. 
The Frcdholm determinants of the even and odd parts 

K?(x, V ) = *^±*^ (272) 
ir(x — y) n(x + y) 

on I/ 2 (0, s; dx) have been studied by Dyson [16]. He used the second derivatives of their logarithms 
to construct potentials for Schrodingcr equations on the half-line, and studied their asymptotics 
with the tools of scattering theory. Jimbo, Miwa, Mori, and Sato [17] related det(l — K s ) to a 
Painlcvc equation. Widom [34] obtained the leading asymptotics using the Krein continuous analog 
of orthogonal polynomials. Deift, Its, and Zhou [11] justified the Dyson asymptotic expansions 
using tools developed for Riemann-Hilbert problems. Tracy and Widom [32] established partial 
differential equations for the Frcdholm determinants of integral operators arising in the study of the 
scaling limit of the distribution functions of eigenvalues of random matrices. We refer the reader 
to the cited references and to [12] for recent results and we apologize for not providing any more 
detailed information here. 

We have, in the present paper, been talking a lot of scattering and determinants and one might 
wonder if this is not a re-wording of known things. In fact, our work is with the multiplicative 
kernels k(xy), and (direct) reduction to additive kernels would lead to (somewhat strange) git + u) 
kernels on semi-infinite intervals (— oo, log(a)]. So wc arc indeed doing something different; one may 
also point out that the entire functions arising in the present study are not of finite exponential 
type; and the scattering matrices do not at all tend to 1 as the frequency goes to infinity. In the 
case of the cosine and sine kernels the flow of information will presumably go from the additive to 
the multiplicative, as the additive situation is more flexible, and has stimulated the development of 
powerful tools, with relation to Painlevc equations, Riemann-Hilbert problems, Intcgrable systems 
[11]. 
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Nevertheless, one may ask if the framework of reproducing kernels in Hilbert spaces of entire 
functions also may be used in the additive situation. This is the case indeed and it is very much 
connected to the method of Krein in inverse scattering theory, and his continuous analog of or- 
thogonal polynomials (used by Widom in the context of the Dirichlct kernel in [34].) The Gaudin 
identities for convolution kernels ([22, App. A16]) play a role very analogous to the identities in the 
present paper (132a), (132b) involved in the study of multiplicative kernels. Widom in his proof [34] 
of the main term of the asymptotics as s — > +oo studied the Krein functions associated with the 
complement of the interval (—1, +1) and he mentioned the interest of extremal properties. In this 
appendix, I shall point out that the resolvent of the Dirichlet kernel indeed does have an extremal 
property: it coincides exactly (up to complex conjugation in one variable) with the reproducing 
kernel of a certain (interesting) Hilbert space of entire functions. This could be a new observation, 
obviously closely related to the method of Widom [34] . 

The space mPW s we shall use is, as a set, the Paley- Wiener space PW S , but the norm is different: 

mPW s = {f(z) entire of exponential type at most s with ||/|| < oo} 

i/wr* (273) 

-i,D 



Let X 8 (z, w) be the element of mPW s which is the evaluator at z: V/ <E mPW s (/, X s (z, •)) = f(z). 
We shall compare X s (z, w) with the resolvent of the kernel 

D s ( X ,y) = S ^f X -f (274) 
tt(x - y) 

on £ 2 (-l, 1; da;). 

Let / G mPW s . It belongs to PW S so 

,<„ . / /W ^'-'T" * - / (275) 
Jn 

On the other hand: 



2Tri(t — z) J R n(t — z) 



f{z)= (J + J J f(t)X a (z,t)dt (276) 

As /(I) = (jr^ + W) X^Tjdt = (JZ^ + W)X s (z, t) dt one has Y^J) = X s (z, t) 
for t G M. We have for y\ and y% real 



X s { yi ,y 2 )= / X s (y u t)X s (y 2 ,t)dt= / X s {y u t)X s (y 2 , t) dt = X s (y 2 , Vl ) (277) 

i»\(-l,l) JR\(-1,1) 

so more generally A s (zT, z 2 ) = A s (z2, z\). 

We apply (276) to f(z) = ^gl"^ for some y G C: 

sin(s(z — y)) f sin(s(t — y)) 



X s (z,t)dt (278) 

-n{z-y) M(-i,i) nt-y) 
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We apply (275) to f(y) = X s {z,y) for some zeC: 



X s (z,y) = f X s (z,t) 



sin(s(t - y)) 
n(t - y) 



dt 



(279) 



Combining we obtain: 



X s (z,y) 



sin(s(z — 



sm(s(t - y)) 
n(t - y) 



dt 



(280) 



tt(z - y) 



Restricting to y € (—1,1), z = x 6 (—1,1), this says exactly: 



x s{x,y) = R s {x,y) 



(281) 



where R s (x, y) is the kernel of the resolvent: 1 + R s = (1 — D s ) 1 , R s — D s — R S D S . The resolvent 
R s (x, y) is entire in (x, y) and the general formula is thus: 
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